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1. Introduction 

The origin of the concept of magnetic charge, if not the name, goes back to antiquity. 
Certain stones in Magnesia, in Anatoha (Asia Minor), were found to exhibit an attractive 
force on iron particles, and thus was discovered magnetism. (Actually there is an ancient 
confusion about the origin of the name, for it may refer to Magnesia, a prefecture in 
Thessaly, Greece, from whence came the settlers ( "Magnets" ) of the city (or the ruins of 
another city) which is now in Turkey. For a recent discussion on the etymology see pp.) 
Electricity likewise was apparent to the ancients, but without any evident connection to 
magnetism. Franklin eventually posited that there were two kinds of electricity, positive 
and negative poles or charge; there were likewise two types of magnetism, north and 
south poles, but experience showed that those poles were necessarily always associated 
in pairs. Cutting a magnet, a dipole, in two did not isolate a single pole, but resulted in 
two dipoles with parallel orientation; the north and south poles so created were bound 
to the opposite poles already existing P]. This eventually was formalized in Ampere's 
hypothesis (1820): Magnetism has its source in the motion of electric charge. That 
is, there are no intrinsic magnetic poles, but rather magnetic dipoles are created by 
circulating electrical currents, macroscopically or at the atomic level. 

Evidently, the latter realization built upon the emerging recognition of the 
connection between electricity and magnetism. Some notable landmarks along the way 
were Oersted's discovery (1819) that an electrical current produced magnetic forces in 
its vicinity; Faraday's visualization of lines of force as a physical picture of electric 
and magnetic fields; his discovery that a changing magnetic field produces a electric 
field (Faraday's law of magnetic induction, 1831); and Maxwell's crowning achievement 
in recognizing that a changing electric field must produce a magnetic field, which 
permitted him to write his equations describing electromagnetism (1873). The latter 
accomplishment, built on the work of many others, was the most important development 
in the 19th Century. It is most remarkable that Maxwell's equations, written down in 
a less than succinct form in 1873, have withstood the revolutions of the 20th Century, 
relativity and quantum mechanics, and they still hold forth unchanged as the governing 
field equations of quantum electrodynamics, by far the most successful physical theory 
ever discovered. 

The symmetry of Maxwell's equations was spoiled, however, by the absence of 
magnetic charge, and it was obvious to many, including Poincare i3j and Thomson 
jU |3] , the discoverer of the electron, that the concept of magnetic charge had utihty, 
and its introduction into the theory results in significant simplifications. (Faraday jB] 
had already demonstrated the heuristic value of magnetic charge.) But at that time, 
the consensus was clearly that magnetic charge had no independent reality, and its 
introduction into the theory was for computational convenience only [Jj, although Pierre 
Curie [8J did suggest that free magnetic poles might exist. It was only well after the 
birth of quantum mechanics that a serious proposal was made by Dirac jH] that particles 
carrying magnetic charge, or magnetic monopoles, should exist. This was based on his 
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observation that the phase unobservabihty in quantum mechanics permits singularities 
manifested as sources of magnetic fields, just as point electric monopoles are sources of 
electric fields. This was only possible if the product of electric and magnetic charges 
was quantized. This prediction was an example of what Gell-Mann would later call the 
"totalitarian principle" - that anything which is not forbidden is compulsory ■ Dirac 
eventually became disillusioned with the lack of experimental evidence for magnetic 
charge, but Schwinger, who became enamored of the subject around 1965, never gave 
up hope. This, in spite of his failure to construct a computationally useful field theory 
of magnetically charged monopoles, or more generally particles carrying both electric 
and magnetic charge, which he dubbed dyons (for his musing on the naming of such 
hypothetical particles, see [H]). Schwinger's failure to construct a manifestly consistent 
theory caused many, including Sidney Coleman, to suspect that magnetic charge could 
not exist. The subject of magnetic charge really took off with the discovery of extended 
classical monopole solutions of non-Abelian gauge theories by Wu and Yang, 't Hooft, 
Polyakov, Nambu, and others [121 CHI Cll CSl CEl CZI ■ With the advent of grand unified 
theories, this implied that monopoles should have been produced in the early universe, 
and therefore should be present in cosmic rays. (The history of magnetic monopoles 
up to 1990 is succinctly summarized with extensive references in the Resource Letter of 
Goldhaber and Trower jjlSij.) 

So starting in the late 1960s there was a burst of activity both in trying to develop 
the theory of magnetically charged particles and in attempting to find their signature 
either in the laboratory or in the cosmos. As we will detail, the former development was 
only partially successful, while no evidence at all of magnetic monopoles has survived. 
Nevertheless, the last few years, with many years of running of the Tevatron, and on the 
eve of the opening of the LHC, have witnessed new interest in the subject, and new limits 
on monopole masses have emerged. However, the mass ranges where monopoles might 
most likely be found are yet well beyond the reach of earth-bound laboratories, while 
cosmological limits depend on monopole fluxes, which are subject to large uncertainties. 
It is the purpose of this review to summarize the state of knowledge at the present 
moment on the subject of magnetic charge, with the hope of focusing attention on the 
unsettled issues with the aim of laying the groundwork for the eventual discovery of this 
exciting new state of matter. 

A word about my own interest in this subject. I was a student of Julian Schwinger, 
and co-authored an important paper on the subject with him in the 1970s Many 
years later my colleague in Oklahoma, George Kalbfieisch, asked me to join him in a new 
experiment to set limits on monopole masses based on Fermilab experiments |2(H I21j. 
His interest grew out of that of his mentor Luis Alvarez, who had set one of the best 
earlier limits on low- mass monopoles [22 I2ni I2H 1211 1211 • Thus, I believe I possess the 
bona fides to present this review. 

Finally, I offer a guide to the reading of this review. Since the issues are technical, 
encompassing both theory and experiment, not all parts of this review will be equally 
interesting or relevant to all readers. I have organized the review so that the main 
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material is contained in sections and subsections, while the third level, subsubsections, 
contains material which is more technical, and may be omitted without loss of continuity 
at a first reading. Thus in section El sections IH.l.lfnmH constitute a detailed proof of 
the quantization condition, while section 13.21 describes the quantum mechanical cross 
section. 

In this review we use Gaussian units, so, for example, the fine-structure constant is 
a = /he. We will usually, particularly in field theoretic contexts, choose natural units 
where h = c = 1. 



2. Classical theory 

2.1. Dual symmetry 

The most obvious virtue of introducing magnetic charge is the symmetry thereby 
imparted to Maxwell's equations in vacuum, 

V ■ E = Airpe, V ■ B = ATipm, 

VxB = -|-E+— je, -V X E= i|-B + — j„. (2.1) 

cot c cot c 

Here pe, je are the electric charge and current densities, and pm, jm are the magnetic 
charge and current densities, respectively. These equations are invariant under a global 
duality transformation. If £ denotes any electric quantity, such as E, pe, or je, while 
A4 denotes any magnetic quantity, such as B, pm, or j^, the dual Maxwell equations 
are invariant under 

£ M, -£, (2.2a) 

or more generally 

^ ^ ^cos^ + TWsin^, M ^ McosO - £sme, {2.2b) 

where ^ is a constant. 

Exploitation of this dual symmetry is useful in practical calculations, even if there is 
no such thing as magnetic charge. For example, its appearance may be used to facilitate 
an elementary derivation of the laws of energy and momentum conservation in classical 
electrodynamics A more elaborate example is the use of fictitious magnetic currents 
in calculate diffraction from apertures (25. (See also [29].) 



2.2. Angular momentum 



J. J. Thomson observed in 1904 |31 EOl |^ the remarkable fact that a static system 

of an electric (e) and a magnetic (g) charge separated by a distance R possesses an 

angular momentum, see figure H The angular momentum is obtained by integrating 

the moment of the momentum density of the static fields: 

\a ^ E X B 
(dr) r X 



Anc 



(dr) r X G 
(dr) r X 



Aire 



er q(r — R) 

^ X ^ 



R 



c 



(2.3) 
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R 

Figure 1. Static configuration of an electric change and a magnetic monopole. 

which follows from symmetry (the integral can only supply a numerical factor, which 
turns out to be in [27J). The quantization of charge follows by applying semiclassical 
quantization of angular momentum: 

J-R= — = n^, n = 0, ±1, ±2, (2.4a) 

or 

eg = m'hc, m' = —. (2.46) 

(Here, and in the following, we use m' to designate this "magnetic quantum number." 
The prime will serve to distinguish this quantity from an orbital angular momentum 
quantum number, or even from a particle mass.) 

2.3. Classical scattering 

Actually, earlier in 1896, Poincare j3] investigated the motion of an electron in the 
presence of a magnetic pole. This was inspired by a slightly earlier report of anomalous 
motion of cathode rays in the presence of a magnetized needle [321 • Let us generalize 
the analysis to two dyons (a term coined by Schwinger in 1969 [TT]) with charges ei, gi, 
and 62, g2, respectively. There are two charge combinations 

q = eie2 + gig2, = • (2.5) 

c 

Then the classical equation of relative motion is [fi is the reduced mass and v is the 
relative velocity) 

d^ r r , , 

/i— r = q— - Kv X —. (2.6) 
at r'^ r-^ 

The constants of the motion are the energy and the angular momentum, 

1 Q 

E = -flV^ + -, J = T X fl\ + KT. (2.7) 

2 r 

Note that Thomson's angular momentum (j2.3|) is prefigured here. 

Because J ■ f = k, the motion is confined to a cone, as shown in figure |21 Here the 
angle of the cone is given by 

X I 

77 = TT' ^ = A^^o^, (2.8) 
Z \k\ 

where vq is the relative speed at infinity, and h is the impact parameter. The scattering 
angle 9 is given by 



d V 



f/2 




Figure 2. The relative motion of two dyons is confined to tlie surface of a cone about 
the direction of the angular momentum. 




Figure 3. Scattering angle as a function of the impact parameter variable x- Here 
kiT = \n\vo/q. 



where 



^ J arctan cot | j , g > 0, 
2 \ ■''^ ~ arctan {^-^^ cot | j , g < 0. 



(2.96) 



When g = (monopole-electron scattering), ^ = vr. The impact parameter h{6) is a 
multiple- valued function of ^, as illustrated in figure El The differential cross section is 
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Figure 4. Classical cross section for monopole-electron and dyon-dyon scattering. 
Again, kir = \K\vQ/q, while g{9) — {^wq / k)'^ {da / dil) . 



therefore 



da 



bdb 



d(cos ( 



— ) E 



4sin^f 



sin X dx 



sin Q dB 



(2.10) 



Representative results are given in |19|, and reproduced here in figure 01 
The cross section becomes infinite in two circumstances; first, when 

sin^ = (sinx^O), = vr, 

we have what is called a glory. For monopole-electron scattering this occurs for 



Xg 



1.047, 1.318, 1.403, 



The other case in which the cross section diverges is when 



(2.11) 



(2.12) 



(2.13) 
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This is called a rainbow. For monopole-electron scattering this occurs at 

9r = 140.1°, 156.7°, 163.5°, .... (2.14) 
For small scattering angles we have the generalization of the Rutherford formula 

Note that for electron-monopole scattering, ei = e, 62 = 0, gi = 0, g2 = g, this 
cross section differs from the Rutherford one for electron-electron scattering by the 
replacement 

^ ^ ^. (2.16) 

V c 

This is a universal feature which we (and others) used in our experimental analyses. 
3. Quantum theory 

Dirac showed in 1931 [H] that quantum mechanics was consistent with the existence of 
magnetic monopoles provided the quantization condition holds, 

eg = m'hc, (3.1) 

where m' is an integer or an integer plus 1/2, which explains the quantization of electric 
charge. This was generalized by Schwinger to dyons: 

em - e2gi = -m'hc. (3.2) 

(Schwinger sometimes argued that m' was an integer, or perhaps an even integer.) 
We will demonstrate these quantization conditions in the following. Henceforth in this 
section we shall set h = c = 1. 

3.1. Vector potential 

One can see where charge quantization comes from by considering quantum mechanical 
scattering. To define the Hamiltonian, one must introduce a vector potential, which 
must be singular because 

V-B^O^B^VxA. (3.3) 

For example, a potential singular along the entire line n is 

gl/nxr nxr\ q ^ 
A(r) = --- z : = -- cot 9(f) (3.4) 



r 2 \r — n ■ r r + n- ry r 
where the latter form applies if n = z, which corresponds to the magnetic field produced 
by a magnetic monopole at the origin, 

Bir) = g^. (3.5) 

In view of fl3.3p . we can write 

B(r) = V X A(r) + ^f(r), (3.6a) 
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where 

V-f(r) =47r5(r), (3.65) 

in which f has support only along the line n passing through the origin. The line of 
singularities is called the string, and f is called the string function. Invariance of the 
theory (wavefunctions must be single-valued) under string rotations implies the charge 
quantization condition 1)3.111 . This is a nonperturbative statement, which is proved in 
section 13.1.31 

3.1.1. Yang's approach Yang offered another approach, which is fundamentally 
equivalent [SH ESI IHEl IHZl IHHl IHOl UHl HI] • He insisted that there be no singularities, 
but rather different potentials in different but overlapping regions: 

= (l-cos^) = ^tan-, 0<7r, (3.7a) 

^ rsin^^ ' r 2' ' ^ ' 

A\ = ^(l + cos^) = --cot-, ^>0. (3.75) 

^ r sin d r 2 

These correspond to the same magnetic field, so must differ by a gradient: 

K-K = Z^n^ = d,\ (3.8) 



where A = 2g(f). Requiring now that e'^'^ be single valued leads to the quantization 
condition, eg = m', m' a half integer. 

3.1.2. Spin approach There is also a intrinsic spin formulation, pioneered by Goldhaber 
nail. The energy 

E = + -, q = 6162 + gm, (3.9) 

2 r 

differs by a gauge transformation from 

„ . 1 U , Zl^^ ^ £. (3.10, 



where 



J = rxp + S, (3.11a) 

/iv = p + ^^. (3.116) 

The quantization condition appears as 

S-r = m'. (3.12) 

The elaboration of this |13] is given in section 13.1.51 
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3.1.3. Strings Let us now discuss in detail the nonrelativistic, quantum scattering 
of two dyons, with electric and magnetic charges ei,gi and e2,g2, respectively. The 
Hamiltonian for the system is 

'H ^ Iruiivf + l-m2vl + - — r, (3.13) 

2 2 ~ ^2] 

where, in terms of the canonical momenta, the velocities are given by 

miVi = pi - eiAe2(ri, t) - giAm2iri, t), (3.14a) 
m2V2 = P2 - e2Aei(r2, t) - 5'2A^i(r2, t). (3.146) 

The electric (e) and magnetic (m) vector potentials are 

47rAe(r, t) = 47rVAe(r, t) - J (dr')f (r - r') x B(r', t), (3.15a) 



47rA„(r, t) = 47rVA^(r, t) + J (dr')*f (r - r') x E(r', t), (3.156) 



with 



Ae(r, t)^ j (dr')f (r - r') • Ae(r', t), (3.16a) 

A„(r,t) = j (dr')*f(r - r') • A^(r',t). (3.166) 

Here, the functions f and *f represent the strings and must satisfy 

V- Wf(r-r') =47r(5(r-r'). (3.17) 

A priori, f and *f need not be related, and could be different for each source. So, for the 
case of dyon-dyon scattering, it would seem that four independent strings are possible. 
The first condition we impose on the Schrodinger equation, 

= E^, (3.18) 

is that it separates when center-of-mass and relative coordinates are employed, which 
implies 

eiAe2(ri,t) = -g2Ami{r2,t) = eig2A{r), (3.19a) 
e2Aei(r2, t) = -^iA^2(ri, t) = e2giA{r), (3.196) 

where r = ri — r2. Correspondingly, there are relations between the various string 
functions, 

*fi(x) = -f2(-x), *f2(x) = -fi(-x), (3.20) 

leaving only two independent ones. The Hamiltonian for the relative coordinates now 
reads 

'H = ^[p- e,g2A{r) + e2giA{r)]^ + ^, (3.21) 
2/1 r 

where /i is the reduced mass. 
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If we furtlier require tliat only one vector potential be present, A. = A', so that 
only the antisymmetric combination of electric and magnetic charges occurring in ()3.2p 
appears, one more relation is obtained between the two f functions, 

f2(x) = -fi(-x). (3.22) 

Notice that ()3.22|1 possesses two types of solutions. 

• There is a single string, necessarily infinite, satisfying 

f(x) = -f(-x). (3.23) 

Then it is easily seen that the vector potential transforms the same way as 
charges and currents do under duality transformations ()2.2&|) . This is the so-called 
symmetric case. 

• There are two strings, necessarily semi-infinite, which are negative reflections of 
each other. 

If identical semi-infinite strings are employed, so that A. ^ A!, the individual charge 
products eig2 and 62(71 occur in the dynamics. The singularities of A. and A! lie on lines 
parallel and antiparallel to the strings, respectively. We will see the consequences for 
the charge quantization condition of these different choices in the following. 

For now, we return to the general situation embodied in ()3.2H1 . For simplicity, we 
choose the string associated with A to he a straight line lying along the direction n, 

1 n X r 

— , semi-mnmte 

r r — (fi ■ r) , , 

11/-^ . . . 3.24 

11/nxr nxr\ 

innnite 



r 2 \r — (n ■ r) r + (fi ■ r 

This result is valid in the gauge in which Xe{m) [ ()3.1(iaj) . ()3.166j) ] is equal to zero. Without 
loss of generality, we will take A! to be given by ()3.24|) with fi — z, which corresponds 
to taking the string associated with A! to point along the —z axis, fi cx — z. 

We now wish to convert the resulting Hamiltonian, 7i, into a form Ti' in which all 
the singularities lie along the z axis. It was in that case that the Schrodinger equation 
was solved in [inj, as described in section 13. 2[ yielding the quantization condition 
()3.2|) . This conversion is effected by a unitary transformation (essentially a gauge 
transformation) , 

n' = e'^He-'^. (3.25) 
The differential equation determining A is 

VA = ei(72[A'(r)- A(r)]. (3.26) 
We take fi to be given by 

fi = sin X cos ipH + sin x sin ipy + cos xz, (3.27) 
and use spherical coordinates [r = (r, 9, (f))], to find 

A= -ei(/2/3(n,r), (3.28) 
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where, for the semi-infinite string (Dirac) 

Pn = (p-i^ + (cos 9 - cos x)F-iO, - V^, x) - 27rr/(x - v), (3.29a) 
and for the infinite string (Schwinger) 

Ps = ^ [(cos e - COS x)F- {6, (P-^,x) + (cos 6 + cos x)F+{e, - V^, x) - 2n7]{x " ^)] • 



(3.296) 



The functions occurring here are 

d(f)' 



F±ie,a,x)-- 
2e(a] 



1 ± cos X cos ^ ± sin X sin 9 cos 0' 

1 ± cos(x + ^^) \ """^^ |q;| 



cos 9 ± cos X 



arctan 



tan ■ 



(3.30) 



l±cos(x-^)y 2 

where the arctangent is not defined on the principal branch, but is chosen such that 
F±{9,a,x) is a monotone increasing function of a. The step functions occurring here 
are defined by 

viO = i'' ^ ^ °' (3.31a) 




^ > 0, 

e(0 = <^ (3.31&) 

The phases, Pd and Ps, satisfy the appropriate differential equation ()3.26|) . for 6* 7^ x 
(as well as 6* 7^ TT — X for Ps) and are determined up to constants. The step functions 
rj are introduced here in order to make e'^ continuous at 6' = x and vr — x, as will be 
explained below. We now observe that 

2tt 

F±{9, 27r + a, x) - F±{9, a, x) = 1 ^- r, (3.32) 

I cost^ ± cosxl 

so that the wavefunction 

^ = e"'^^', (3.33) 

where is the solution to the Schrodinger equation with the singularity on the z axis, 
is single-valued under the substitution — >■ + 27r when the quantization condition 
()3.1|) is satisfied. 

Notice that integer quantization follows when an infinite string is used while a 
semi-infinite string leads to half-integer quantization, since Ps changes by a multiple 
of 271 when — + 27r, while Pd changes by an integer multiple of An. Notice that 
Pd possesses a discontinuity, which is a multiple of 47r, at 6* = x, while Ps possesses 
discontinuities, which are multiples of 27r, at 9 = x^ ~ X- virtue of the above- 
derived quantization conditions, e'^ is continuous everywhere. Correspondingly, the 
unitary operator e^^, which relates solutions of Schrodinger equations with different 
vector potentials, is alternatively viewed as a gauge transformation relating physically 
equivalent descriptions of the same system, since it converts one string into another. 
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[Identical arguments applied to the case when only one vector potential is present leads 
to the condition ()3.2j) . where m' is an integer, or an integer plus one-half, for infinite 
and semi-infinite strings, respectively.] 

It is now a simple application of the above results to transform a system 
characterized by a single vector potential with an infinite string along the direction 
n into one in which the singularity line is semi-infinite and lies along the +z axis. This 
can be done in a variety of ways; particularly easy is to break the string at the origin 
and transform the singularities to the z axis. Making use of ()3.28j) with eig2 — > — m'/2 
and ()3.29aj) for n and — n, we find 

A = m'P's{n,r) with p'^ = (J) - ip + p^. (3.34) 

In particular, we can relate the wavefunctions for infinite and semi-infinite singularity 
lines on the z axis by setting x = in ()3.34|1 . 

(3's = (P-^, (3.35) 

so 

^(infinite) = e-'"''^'^-'^)^ (semi-infinite). (3.36) 
Note that ()3.36|1 or ()3.34j) reiterates that an infinite string requires integer quantization. 

3.1.4- Scattering In the above subsection, we related the wavefunction when the string 
lies along the direction n with that when the string lies along the z axis. When there 
is only a single vector potential (which, for simplicity, we will assume throughout the 
following), this relation is 

where /3 is given by ()3.29a|l . (j3.29&j) . or ()3.34|1 for the various cases. For concreteness, 
if we take to be a state corresponding to a semi-infinite singularity line along the 
+z axis, then /? is either /Sd ()3.29a|l or P'g ()3.34|1 depending on whether the singularity 
characterized by n is semi- infinite or infinite. By means of ()3.37j) . we can easily build 
up the relation between solutions corresponding to two arbitrarily oriented strings, with 
n and n' say, 

which expresses the gauge covariance properties of the wavefunctions. 

For scattering, we require a solution that consists of an incoming plane wave and 
an outgoing spherical wave. We will consider an eigenstate of J • k where J is the total 
angular momentum 

J = r X (p + m'An) + m'r, (3.39) 

and k is the unit vector in the direction of propagation of the incoming wave (not 
necessarily the z axis). This state cannot be an eigenstate of k- (r x p), since this 
operator does not commute with the Hamiltonian. However, since 

e'^k ■ Je"'^ = k ■ (r X p) - m', (3.40a) 
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for a reorientation of the string from n to k, 

A = m'[/3(n,r)-/3B(k,r)], (3.406) 

because 

k ■ (f X (k X f)) - , , 

^ ^ ^ = l + k-f, (3.41) 

1 — k f 

the incoming state with eigenvalue ^H] 

(k ■ J)' = -m' (3.42) 
is simply related to an ordinary modified plane wave [77 is defined below in fl3.47|) ] 

^.^ = e"'^ exp {i [k ■ r + 7/ ln(A;r - k ■ r)]} . (3.43) 

This state exhibits the proper gauge covariance under reorientation of the string. 
The asymptotic form of the wavefunction is 

^ _ g-W/3(fi,r) A^._^yrn^(^^ym'^^ ^-^ - r]\n2kr --L + 6l), r cx).(3.44) 

jm 

The summation in ()3.44|) is the general from of the solution when the singularity line is 
semi-infinite, extending along the +z axis. In particular, is a generalized spherical 
harmonic, which is another name for the rotation matrices in quantum mechanics 

{jm'\e'^'^^e''^-^^e'^^^\jm) = e'"^'^^^=V™;(f) = e''^'^U^JJ(e)e''^^ , (3.45) 

V2j + 1 

6l is the Coulomb phase shift for noninteger L, 

5l = argr(L+ l + iT]), (3.46) 

and 



^ + 1 = \l y + 1) - m'^, ^ = y ' ^ = ^1^2 + 9i92. (3.47) 
Upon defining the outgoing wave by 

^ ~ e"^^ (gi[k.r+r,ln(fcr-k.r)] ^ ^^^^^ ^ ^g ^g^, 

where A is given by ()3.40 we find that 

*out = le'('="-''''^'='-)e^™'^/(^). (3.49) 
r 

In terms of the scattering angle, 6, which is the angle between k and r, the scattering 
amplitude is 

00 

2iA;/(^") = ^2^iy%,,{'^ - 9, 0)e-'(-^-2^-). (3.50) 

j=\m'\ 

The extra phase in ()3.49j) is given by (where k characterized by 9' , 0' and — k by tt — 0', 

0'±7r) 

7 = /3D(k, -k) + - 0' T TT - /3z5(k, r) + 0, (3.51) 
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where 



^ _ cos (^±^) sin ' 



arctan -6 = — — -^7- , , ; . (3.52) 

2^ cos (^^) cos (^^) ^ ^ 

Straightforward evaluation shows that 

I = (mod27r), (3.53) 
so that there is no additional phase factor in the outgoing wave. 

3.1.5. Spin Classically, the electromagnetic field due to two dyons at rest carries 
angular momentum, as in section 12. 2| 

Sclassical = ^'y. (3.54) 

A quantum-mechanical transcription of this fact allows us to replace the nonrelativistic 
description explored above, in which the interaction is through the vector potentials 
(apart from the Coulomb term), by one in which the particles interact with an intrinsic 
spin. The derivation of the magnetic charge problem from this point of view seems first 
to have carried out by Goldhaber |321 in a simplified context, and was revived in the 
context of 't Hooft-Polyakov monopoles [TSl CH HZl EHl 113 IIHI , where the spin is called 
"isospin." 

Before introducing the notion of spin, we first consider the angular momentum of the 
actual dyon problem. For simplicity we will describe the interaction between two dyons 
in terms of a single vector potential A., and an infinite string satisfying ()3.23p . (The 
other cases are simple variations on what we do here, and the consequences for charge 
quantization are the same as found in section Til. 31 ) Then the relative momentum of 
the system is 

p = /iv — m'A.. (3.55) 

Since from ()3.15ajl . ()3.15&|1 . we have the result in ()3.(ia|l . or 

Vx^=^-f(r), (3.56) 

we have the following commutation property valid everywhere, 

r 



yuv X yuv = — im' 



ffr) 



(3.57) 

Motivated by the classical situation, we assert that the total angular momentum 
operator is ()3.39|) . or 

J = r X /iv + m'i. (3.58) 

This is confirmed ^T] by noting that, almost everywhere, J is the generator of rotations: 

y[r, J ■ H = •^'^ X r, (3.59a) 

T[/iv, J ■ 5u] = 6uj X — m'f (r) x {6u} x r), (3.596) 
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where 6ljJ stands for an infinitesimal rotation. The presence of the extra term in ()3.595|) 
is consistent only because of the quantization condition [32] • For example, consider the 
effect of a rotation on the time evolution operator, 



e-'-^-^'^exp 



-i / dtn 



A3-5u} 



exp 



— 1 



where 



dt {n + 5n) 



5y = 5ijj ^ r. 



5n = i[l-LJ- 5u:] = m'v ■ [f (r) x 5r], 
Using the representation for the string function, 

f (r) = 47r /" dxi [(5(r - x) - 6{r + x)] 
Jc 2 

where C is any contour starting at the origin and extending to infinity, and the notation 
dtv = dr, we have 



(3.60) 
(3.61) 
(3.62) 



i j dtSH = -im'An j dr ■ (dx x 5r)^ [5(r - x) - 5{y + x)] . 



(3.63) 



Since the possible values of the integral are 0, ±|, ±1, the unitary time development 
operator is unaltered by a rotation only if m' is an integer. (Evidently, half-integer 
quantization results from the use of a semi-infinite string.) 

Effectively, then, J satisfies the canonical angular momentum commutation 
relations (see also section in.l.6j) 
1 



-J X J = J, 
1 

and is a constant of the motion, 



0. 



(3.64) 



(3.65) 



And, corresponding to the classical field angular momentum ()3.54j) . the component of 
J along the line connecting the two dyons, m', should be an integer. 

The identification of m' as an angular momentum component invites us to introduce 
an independent spin operator S. We do this by first writing [TT], as anticipated in ()3.12|) . 

m' = S ■ f , (3.66a) 

and 

S X r 



/iv = p + 



(3.666) 



which, when substituted into ()3.58|1 yields ()3.11ajl . 

J = rxp + S. (3.67) 

We now ascribe independent canonical commutation relations to S, and regard ()3.66aj) 
as an eigenvalue statement. The consistency of this assignment is verified by noting 
that the commutation property 



/iv X /iv = — im 



(3.68) 
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holds true, and that S ■ f is a constant of the motion, 
[S-f,/iv] = 0. 

In this angular momentum description, the Hamiltonian, 
form 



n 



1 

2/i 



p" + 



2S ■ L 



+ 



fS-f) 



+ 



in terms of the orbital angular momentum, 
L = r X p. 

The total angular momentum J appears when the operator 



(rp) +Yrp 



is introduced into the Hamiltonian 



n = — 

2/i 



Pr 



iJ-r) 



1 

r 
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(3.69) 

■211) . can be written in the 
(3.70) 

(3.71) 

(3.72) 

(3.73) 



fJ-f) 



m 



(3.74) 



In an eigenstate of and J ■ f . 

(J7=J(J + 1), 

()3.73|) yields the radial Schrodinger equation ()3.128aj) solved in section 13.21 This 
modified formulation, only formally equivalent to our starting point, makes no reference 
to a vector potential or string. 

We now proceed to diagonalize the S dependence of the Hamiltonian, ()3.7()j) or 
()3.73|) . subject to the eigenvalue constraint 

(S • f )' = m'. (3.75) 

This is most easily done by diagonalizing [IH] the angular momentum operator ()3.(i7|l . In 
order to operate in a framework sufficiently general to include our original symmetrical 
starting point, we first write S as the sum of two independent spins 

S = Sa + S5. (3.76) 

We then subject J to a suitable unitary transformation j32] 

3' = U3U-\ (3.77) 

where 

U = exp[i(S, ■ 4>)0] exp[i(Sfe ■ 0)(0 - vr)], (3.78) 

which rotates Sa,b ■ ? into ±(Sa,fe)3. This transformation is easily carried out by making 
use of the representation in terms of Euler angles, 

exp(iS • 4>6) = exp(— 105*3) exp(i6'S'2) exp(i0S'3). (3.79) 

The general form of the transformed angular momentum, 

Sal, , 



J' 



r X 



p H sm u 

r 



1 + cos 6* 



+ 



cos 6' 



+ f(S„-S„)J 



3' 



(3.80) 
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is subject, a priori, only to the constraint (j3.75|) . or 

(S, - S,y, = m'. (3.81) 

We recover the unsymmetrical and symmetrical formulations by imposing the following 
supplementary eigenvalue conditions: 

(1) : S',, = 0, (3.82a) 

(2) : (S, + S,y, = 0. (3.826) 

These yield the angular momentum in the form ()3.58p or ()3.39|) . the vector potential 
appearing there being, respectively, 

(1) : A=--^cot-, (3.83a) 

(2) : A=— ^cot^, (3.836) 

which are (HT^ with n = z. [See also (1X71 . (Q-] 

The effect of this transformation on the Hamiltonian is most easily seen from the 
form (IT7H|) . 



U 



Pr + - 



t/-^=P.' + ^(rx/iv)2 = (/xv)^ (3.84) 



making use of ()3.72|1 . or 



n' = UnU~^ = -nv^ + -. (3.85) 
2 r 

So by means of the transformation given in ()3.78|) we have derived the explicit 
magnetic charge problem, expressed in terms of J' and Ti' , from the implicit formulation 
in terms of spin. These transformations are not really gauge transformations, because 
the physical dyon theory is defined only after the eigenvalue conditions 1)3.811) and 
(I3.82a|l - ()3.826jl are imposed. The unsymmetrical condition (1), ()3.82ajl . gives rise to 
the Dirac formulation of magnetic charge, with a semi-infinite singularity line, and, 
from ()3.81)1 . m' either integer or half- integer. The symmetrical condition (2), ()3.826)1 
gives the Schwinger formulation: An infinite singularity line [with ()3.23|) holding], 
and integer quantization of m' . These correlations, which follow directly from the 
commutation properties of angular momentum (the group structure), are precisely the 
conditions required for the consistency of the magnetic charge theory, as we have seen 
in section in. 1.3[ 

Even though the individual unitary operators U are not gauge transformations, a 
sequence of them, which serves to reorient the string direction, is equivalent to such a 
transformation. For example, if we formally set = in ()3.78p . 

?7(i) = exp(iS -0(0-71)], (3.86) 

we have the transformation which generates a vector potential with singularity along 
the positive z axis, ()3.83aj) . while 

[/(2) = exp[iS ■ U2(e - tt)] (3.87) 
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generates a vector potential with singularity along n, the first form in (j3.24p . where 
is the angle between n and r, 



cos G = cos 9 cos X + sin 6^ sin x cos{(f) — ip) 
[the coordinates of n are given by ()3.27|) ]. and 



U2 



n X r 

In X r| 



The transformation which carries (|3.83ap into the first form in (|3.24|) is 



(12) 



(3.88) 
(3.89) 
(3.90) 



Since f/(i2) reorients the string from the direction z to the direction fi, it must have the 
form 



(12) 



exp(iS ■ n<l>) exp(— iS ■ tpx)- 



(3.91) 



The angle of rotation about the n axis, $, is most easily determined by considering the 
spin-1/2 case S = |cr, and introducing a right-handed basis, 



n X r 



Ul 



n, 



U2 



|n X r| 

Then straightforward algebra yields 



U3 = n X U2. 



cos -$ 
2 



sin -$ 
2 



1 sin ^6' cos |x ~ cos |6' sin cos((/) — ^Z') 



sin |9 



cos ^9 sin |x sin(0 — ip) 
sin 1 6 



(3.92) 



(3.93a) 



(3.936) 



The corresponding transformation carrying the vector potential with singularities along 
the negative z axis [ ()3.83ajl with 9 —* 9 — tt], into the vector potential with singularities 
along the direction of — fi [the first form in ()3.24|) with n —* — n], are obtained from 
flCT]) and ^M^-^MB [see also and (jSSZI)] by the substitutions 

9^9 + 7T, 0-^6 + vr. (3.94) 

The combination of these two cases gives the transformation of the infinite string, of 
which (j3.78j) is the prototype. 

Since the effect of exp(— iS ■ ipx) is completely given by 

exp(-iS • MS3 exp(iS ■ V'x) = S ■ fi, (3.95) 

that is, for the transformation ()3.9H1 . 





(12) 



r X I p + cot -S3 1 - rS^ 



^(12) 



exp(iS • n$) 



^ 
r X I p H cot -b ■ n 

r 2 



f S • fi 



exp(-iS • n<l)), (3.96) 
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in a state when S ■ n has a definite eigenvalue — m', [/(12) is effectively just the gauge 
transformation which reorients the string from the z axis to the direction n. And, 
indeed, in this case, 

i$ = i/?z5 (mod27r), (3.97) 

where /3_d is given by ()3.29ajl as determined by the differential equation method. 

3.1.6. Singular gauge transformations We now make the observation that it is precisely 
the singular nature of the gauge transformations ()3.25|) and ()3.85|) which is required for 
the consistency of the theory, that is, the nonobservability of the string. To illustrate 
this, we will consider a simpler context, that of an electron moving in the field of a static 
magnetic charge of strength (7, which produces the magnetic field 

B = g^. (3.98) 

The string appears in the relation of B to the vector potential, ()3.56|) . ()3.6a|) . or 

B = VxA + ^f(r), (3.99) 

where the string function f satisfies ()3.17|) . Reorienting the string consequently changes 
A, 

A ^ A', (3.100) 
which induces a phase change in the wavefunction, 

m ^m' = e'^m. (3.101) 

The equation determining A is (j3.26|) . or 

VA = e(A'-A), (3.102) 
which makes manifest that this is a gauge transformation of a singular type, since 

V X VA ^ 0. (3.103) 

Recognition of this fact is essential in understanding the commutation properties 
of the mechanical momentum (called /iv above), 

7r = p-eA, (3.104) 

since 

77 X TT = -V X V + ie(V X A). (3.105) 

(Here, the parentheses indicate that V acts only on A, and not on anything else to the 
right.) Consider the action of the operator ()3.105p on an energy eigenstate Certainly 
V X = away from the string; on the string, we isolate the singular term by making 
a gauge transformation reorienting the string, 

^ = e"'^^', (3.106) 
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where is regular on the string associated with A. Hence 

{0 off string, 

(3.107) 
i(V X VA)\E' on string, 

so by jSini and 

- V X V^(r) = ie^f(r)^(r). (3.108) 

Thus, when acting on an energy eigenstate [which transforms hke ()3.101|) under a string 
reorientation], ()3.105|1 becomes 

77 X TT ^ ie[(V X A) + ^f(r)] = ieB. (3.109) 

This means that, under these conditions, the commutation properties of the angular 
momentum operator ()3.58j) . 

J = rX7r-e^f, (3.110) 

are precisely the canonical ones 

^[r,3 -Su)]^ 6u) xr, (3.111a) 

-[7r,J ■ 6uj] ^ 6uj X TV. (3.111&) 

In section 13.1.51 we considered the operator properties of J on the class of states for 
which V X V = 0, so an additional string term appears in the commutator ()3.596|1 . 
Nevertheless, in this space, J is consistently recognized as the angular momentum, 
because the time evolution operator is invariant under the rotation generated by J. Here, 
we have considered the complementary space, which includes the energy eigenstates, in 
which case the angular momentum attribution of J is immediate, from ()3.111aj) - ()3.1115|) . 

Incidentally, note that the replacement (j3.109|) is necessary to correctly reduce the 
Dirac equation describing an electron moving in the presence of a static magnetic charge, 

(77r + m)^ = 0, (3.112) 

to nonrelativistic form, since the second-order version of ()3.112|1 is 

(vr^ + m^-eo- -6)^ = 0, (3.113) 

where B is the fully gauge invariant, string independent, field strength ()3.98|1 . rather 
than (V X A), as might be naively anticipated. This form validates the consideration 
of the magnetic dipole moment interaction, including the anomalous magnetic moment 
coupling, which we will consider numerically below, both in connection with scattering 
(section [3.2. ip and binding (section E)). 

Similar remarks apply to the non-Abelian, spin, formulation of the theory, given by 
()3.70|) . If we define the non-Abelian vector potential by 

eA = -^, (3.114) 
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the mechanical momentum ()3.666|) of a point charge moving in this field is again given 
by dSini, or 

7r = p-eA, (3.115) 

and the magnetic field strength is determined, analogously to ()3.109|) . by 

1 f 
eB = -TT X 77 = (V X eA) - ieA x eA = -S • f— . (3.116) 
i 

This reduces to the Abelian field strength ()3.98|) in an eigenstate of S ■ f , 

(S-f)' = -e^, (3.117) 

which is a possible state, since S ■ f is a constant of the motion, 

[S-f,7r]=0. (3.118) 

The Abelian description is recovered from this one by means of the unitary 
transformation ()3.79p 

U = exp{-i(f)S3) exp{ieS2) exp(i053). (3.119) 

Under this transformation, the mechanical momentum, (j3.115|l . takes on the Abelian 
form, 

t/7r[/-i = p + 0— tan-, (3.120) 
r 2 

where we see the appearance of the Abelian potential 

d) 

eA = -53-tan-, (3.121) 
r 2 

corresponding to a string along the —z axis. In an eigenstate of S3, 

S's = {US ■ rU-^y = -eg, (3.122) 

this is the Dirac vector potential ()3.7ajl . To find the relation between this vector potential 
and the field strength, we apply the unitary transformation ()3.120|) to the operator 

eB = V X eA + eA X V-ieA X eA (3.123) 
to obtain, using Stokes' theorem, 

UeBU-^ = (V X eA) - if/V x Vf/"^ = (V x eA) - Ssfir), (3.124) 
where f is the particular string function 

f(r) = -Ankri{-z)S{x)S{y), (3.125) 
r] being the unit step function. In this way the result ()3.99p is recovered. 
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3.1.7. Commentary There is no classical Hamiltonian theory of magnetic charge, since, 
without introducing an arbitrary unit of action [1^11201 , unphysical elements (strings) 
are observable. In the quantum theory, however, there is a unit of action, h, and 
since it is not the action W which is observable, but exp{iW/K), a well-defined theory 
exists provided charge quantization conditions of the form (I3.2|l or ()3.1|1 are satisfied. 
The precise form of the quantization condition depends on the nature of the strings, 
which define the vector potentials. It may be worth noting that the situation which first 
comes to mind, namely, a single vector potential with a single string, implies Schwinger's 
symmetrical formulation with integer quantization jHSj. 

We have seen in the nonrelativistic treatment of the two-dyon system that the 
charge quantization condition is essential for all aspects of the self-consistency of the 
theory. Amongst these we list the nonobservability of the string, the single-valuedness 
and gauge- covariance of the wavef unctions, and the compatibility with the commutation 
relations of angular momentum. In fact, all these properties become evident when it 
is recognized that the theory may be derived from an angular momentum formulation 

3.2. Nonrelativistic Hamiltonian 

We now must turn to explicit solutions of the Schrodinger equation to obtain numerical 
results for cross sections. For a system of two interacting dyons the Hamiltonian 
corresponding to symmetrical string along the entire z axis is 

h"^ 2m' cos 6' 1 d m'^ •? ^\ Q , . 

H = -— V2 + — -— cot^^^ 3.126 
2fi \ sm 6 1 0(f) J r 

where the quantity k in ()2.5p is replaced by the magnetic quantum number m' defined in 
()3.2|) . [This is ()3.2H) with A.= ^ given by the second form in ()3.24|) with fi = z.] Even 
though this is much more complicated than the Coulomb Hamiltonian, the wavefunction 
still may be separated: 

^(r) = R{r)e{9)e'""t', (3.127) 

where the radial and angular factors satisfy 

^^l^^t^.p.M±llzI^)u = 0, (3.128.) 



1 d / ^ d \ m — 2mm' cos 9 + m 
sin 9- 



■n.n\—.nl -2. Q = J {j + 1)Q .{3.128b) 

sm9d9 \ d9 J sin^ 9 

The solution to the 9 equation is the rotation matrix element: [x = cos 9) 

U^llid) = (jm'le'-^^^/^ljm) oc (1 - x)^(l + x)^F,(l-'-'-^-)(x), (3.129) 

where pj™'"^ are the Jacobi polynomials, or "multipole harmonics" [lU]- This forces m' 
to be an integer. The radial solutions are, as with the usual Coulomb problem, confluent 
hypergeometric functions, 

R^.(r) = e-'^''{kr)^F{L + 1 - ir/, 2L + 2, 2iA;r), (3.130a) 
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V=g-^, ^ = ^' L + l = ^l{j + -^] -m-. (3.1306) 

Note that in general L is not an integer. 

We solve the Schrodinger equation such that a distorted incoming plane wave is 
incident, 

^in = exp{i [k ■ r + ?7ln(A;r - k ■ r)]} . (3.131) 
Then the outgoing wave has the form ()3.49j) . (here 6 is the scattering angle) 

r 

where the scattering amplitude is given by ()3.5()j) . or 

00 

2ikf{9) = (2j + l)uilU^ - 9)e--^^-^-''^^ (3.133) 

j=\m'\ 

in terms of the Coulomb phase shift ()3.4(ij) . 

5L = argr(L + l + ir/). (3.134) 

Note that the integer quantization of m' results from the use of an infinite ( "symmetric" ) 
string; an unsymmetric string allows m' = integer + i. 

We reiterate that we have shown that reorienting the string direction gives rise to 
an unobservable phase. Note that this result is completely general: the incident wave 
makes an arbitrary angle with respect to the string direction. Rotation of the string 
direction is a gauge transformation. 

By squaring the scattering amplitude, we can numerically extract the scattering 
cross section. Analytically, it is not hard to see that small angle scattering is still given 
by the Rutherford formula ()2.15|1 : 

for elect ron-monopole scattering. The classical result is good roughly up to the first 
classical rainbow. In general, one must proceed numerically. In terms of 

9iO) = ^,\fm' (3.136) 

we show various results in figures EHZI Structures vaguely reminiscent of classical 
rainbows appear for large m', particularly for negative rj, that is, with Coulomb 
attraction. 

3.2.1. Magnetic dipole interaction We can also include the effect of a magnetic dipole 
moment interaction, by adding a spin term to the Hamiltonian, 

ns = -^^ia.B, B = g^. (3.137) 
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Figure 6. Quantum dyon-dyon scattering, m' = 1. 
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Figure 8. Spinflip (F) and nonflip (NF) cross sections for various gyromagnetic ratios. 
The first graph shows m' = +1, the second m' — —1. 



For small scattering angles, the spin-flip and spin-nonflip cross sections are for 7 = 1 
(0 « 1) 



da 



m 
2k 



sin2^/2 
sin^^/2' 



da 







NF 


\2k) 



cos2 0/2 
sin^ e/2 ' 



(3.138) 



Numerical results are shown in figure |H1 and figure IHl Note from the figures that the spin 
flip amplitude always vanishes in the backward direction; the spin nonflip amplitude also 
vanishes there for conditions almost pertaining to an electron: m' > 0, 7 = 1. 
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Figure 9. Spinflip and nonflip cross sections for m' = 10, 20. 

The calculations shown in figures ISHHl were done many years ago [12], which just 
goes to show that "good work ages more slowly than its creators." The history of the 
subject goes much further back. Tamm [Ml calculated the wavefunction for the electron- 
monopole system immediately following Dirac's suggestion P]; while Banderet |55j . 
following Fierz |56^, was the first to suggest a partial- wave expansion of the scattering 
amplitude for the system. The first numerical work was carried out by Ford and Wheeler 

while the comparison with the classical theory can be found, for example, in jHElEni- 

3.3. Relativistic calculation 

A relativistic calculation of the scattering of a spin-1/2 Dirac particle by a heavy 
monopole was given by Kazama, Yang, and Goldhaber |39J. They used Yang's 
formulation of the vector potential described above in section 13.1.11 In order to arrive 
a result, they had to add an extra infinitesimal magnetic moment term, in order to 
prevent the charged particle from passing through the monopole. The sign of this 
term would have measurable consequences in polarization experiments. It does not, 
however, appear in the differential cross sections. It also does not affect the helicity 
flip and helicity nonflip cross sections which are shown in flgure ^| The vanishing of 
the helicity nonflip cross section in the backward direction precisely corresponds to the 
vanishing of the nonrelativistic spinflip cross section there. The correspondence with 
the nonrelativistic calculation with spin seems quite close. 

4. Non-Abelian monopoles 

Although the rotationally symmetric, static solution of the Yang-Mills equations was 
found by Wu and Yang in 1969 it was only in 1974 when 't Hooft and Polyakov 
included the Higgs field in the theory that a stable monopole solution was found P^IT^. 
Dyonic configurations, that is, ones with both magnetic and arbitrary electric charge. 
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Figure 10. Relativistic helicity-flip and helicity-nonflip cross sections. Note for ( 
helicity nonflip corresponds to spin flip, while helicity flip means spin nonflip. 



were found by Julia and Zee [T^. Here we discuss the unit monopole solution first 
obtained by Prasad and Sommerfield jHOl- Referred to as BPS monopoles, which saturate 
the Bogomolny energy bound [HI], they are static solutions of the SU(2) theory 

^ -FrFa^^u + \{D^HY{D^Hf - ^H'^H'^ - v'')\ (4.1) 



167r "^'^ 2^ ^ ^ ' 4 



where 



D^H = d^H - eA^ X H, (4.2) 

for an isotopic triplet Higgs field H. (We denote the coupling strength by e to avoid 
confusion with the magnetic charge g.) This is the Georgi-Glashow model [HI], in which 
the massive vector boson has mass 



= v47ref, (4.3a) 
while the Higgs boson mass is 

mn = VXv. (4.36) 
This model possesses nontrivial topological sectors. The topological charge is 

k = J (dx)ei,fctr(F,feA^). (4.4) 

In the limit of zero Higgs coupling, A = 0, where the Higgs boson mass vanishes, 
Bogomolny showed that the classical energy of the configuration was bounded by the 
charge. 
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The solution found by Prasad and Sommerfield achieves the lower bound on the energy, 
E = v^ivr/cf/e, with unit charge, k = 1, and has the form 

H = ^-ah{r), eA = o-x^/(r), (4.6a) 



where with = v 47ref r 

%)=^(cothe-0, ^^'") = ^~^' ^^-^^^ 

If we regard (j"/2 as the isotopic generator, the forms of the isotopic components of the 
Higgs field and vector potential are 

H'^ = 2h{r)^, eA^ = 2/(r)e,,,^. (4.7) 

These fields describe a monopole centered on the origin. Far away from that (arbitrary) 
point, the behavior of the fields is given by 

r oo : h{r) v, f{r) 1, (4.8) 

so we see that the vector potential ()4.6(7| indeed describes a monopole of twice the Dirac 
charge, according to ()3.114|1 . because the spin is S = |cr. However, this monopole has 
structure, and is not singular at the origin, because both / and h vanish there. The 
energy of the configuration, the mass of the monopole, is finite, as already noted, 

E = [ (dr)V^ ti H"^ = - [ dS ■ Vh^(r) = Va^—, (4.9) 
J 2 Js^ ve 

which takes into account the relation between r and C,, as claimed. In physical units 
this gives a mass of this monopole solution of 

M=^, (4.10) 
a 

in terms of the "fine structure" constant of the gauge coupling, a = in Gaussian 
units. If the electroweak phase transition produced monopoles, then, we would expect 
them to have a mass of about 10 TeV jHSl lElj- There are serious doubts about this 
possibihty [HSl IHSl EH IHH] so it is far more likely to expect such objects at the GUT 
scale. Magnetic monopole solutions for gauge theories with arbitrary compact simple 
gauge groups, as well as noninteracting multimonopole solutions for those theories have 
been found 

In general, the classical 't Hooft-Polyakov monopole mass in the Georgi-Glashow 
model is with nonzero Higgs mass is 

Mci = fi{z), z = ——. 4.11 

a Mw 

The function /i(0) = 1, and is less than 2 for large z. Quantum corrections to the 
classical mass have been considered on the lattice [70] . 

We will not further discuss non-Abelian monopoles in this review, because that is 
such a vast subject, and there are many excellent reviews such as jZU [Z2], as well as 
textbook discussions [731 • merely note that asymptotically, an isolated non-Abelian 
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monopole looks just like a Dirac one, so that most of the experimental limits apply 
equally well to either point-like or solitonic monopoles. Even the difficulties with the 
Dirac string, which have not been resolved in the second quantized version, to which we 
now turn our attention, persist with non-Abelian monopoles 

5. Quantum field theory 

The quantum field theory of magnetic charge has been developed by many people, 
notably Schwinger [13 EHl EH EHl ES] and Zwanziger [ZOl EH IHD IH2] • We should cite 
the review article by Blagojevic and Senjanovic jHSl; which cites earlier work by those 
authors. A recent formulation suitable for eikonal calculations is given in jHl|, and will 
be described in section 15.31 and following. 

5.1. Lorentz invariance 

Formal Lorentz invariance of the dual quantum electrodynamics system with sources 
consisting of electric charges {ca} and magnetic charges {ga} was demonstrated provided 
the quantization condition holds: 

/ / f , unsymmetric 1 
eaQb - ebga = fn = i ^ }, n E Z. (5.1) 

n , symmetric J 

"Symmetric" and "unsymmetric" refer to the presence or absence of dual symmetry in 
the solutions of Maxwell's equations, refiecting the use of infinite or semi-infinite strings, 
respectively. 

5.2. Quantum action 

The electric and magnetic currents are the sources of the field strength and its dual 
(here, for consistency, we denote by j^, what we earlier called j^, j^, respectively): 

= 47rj^ and 5^*F^, = 47r%, (5.2) 

where 

which imply the dual conservation of electric and magnetic currents, and 
respectively, 

d^j^ = 0, and d^Y = 0. (5.4) 

As we will detail below, the relativistic interaction between an electric and a magnetic 
current is 



W{j, *j) = J (dx)(dx')(da;") T{x)e,uardT {x - x') D+ {x' - x")f {x") . (5.5) 
Here the electric and magnetic currents are 

jf, = eij'y^ifj and % = gxit^X, (5.6) 
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for example, for spin-1/2 particles. The photon propagator is denoted by D^{x — x') 
and /^(x) is the Dirac string function which satisfies the differential equation 

d^r{x)=A'K6{x), (5.7) 

the four-dimensional generalization of ()3.17j) . A formal solution of this equation is given 
by 

/^(x) =4™'^ (n- 9)^^(5 (x), (5.8) 

where is an arbitrary constant vector. [Equation ()3.125|1 results if fi = — z, in which 
case f(r,t) = f(r)5(t).] 

5.3. Field theory of magnetic charge 

In order to facilitate the construction of the dual-QED formalism we recognize that the 
well-known continuous global U(l) (itxa/ symmetry (j2.26jl [ZHl EHl ESj implied by ()5.2|1 . 
()5.4|) . given by 




cos 9 sin 9 

- sin 9 cos 9 

cos 9 sin 9 

— sin 9 cos 9 




(5.9a) 



(5.96) 



suggests the introduction of an auxiliary vector potential B^{x) dual to A^{x). In order 
to satisfy the Maxwell and charge conservation equations, Dirac p3| modified the field 
strength tensor according to 

F^u = d^A, - d,A^ + (5.10) 

where now ()5.2p gives rise to the consistency condition on G^y{x) = —Gy^{x) 

d-'*F^, = -d''G,,=47T*j^. (5.11) 

We then obtain the following inhomogeneous solution to the dual Maxwell's equation 
()5.11|) for the tensor G^,y(x) in terms of the string function and the magnetic current 
%■■ 

G^^{x) = 47r (n ■ 9)~^ [n^*j^{x) - n^*jf,{x)] 

(dy) lUx - y)%iy) - Ux - y)%(y)] , (5.12) 



where use is made of ()5.4|) . ()5.7|) . and ()5.8|) . A minimal generalization of the QED 
Lagrangian including electron-monopole interactions reads 

C = -^F^uF^" + ii-id + e^A -m^)ij + x (^7-9 - x, (5.13) 

where the coupling of the monopole field x{^) to the electromagnetic field occurs 
through the quadratic field strength term according to ()5.1()j) . We now rewrite the 
Lagrangian ()5.13|) to display more clearly that interaction by introducing the auxiliary 
potential B^{x). 
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Variation of (|5.13|) with respect to the field variables, ip, X and A^, yields in addition 
to the Maxwell equations for the field strength, F^j^, ()5.2|) where = e?/'(x)7^^/'(x), 

the equation of motion for the electron field 

(17^ + 67^(0;) - m^)^(a;) = 0, (5.14) 

and the nonlocal equation of motion for the monopole field, 

(i79 - x{x) -l^j (dy) *F^^{y)^^ = . (5.15) 

[We regard Gfj,^{x) as dependent on x, x but not A^. Thus, the dual Maxwell equation 
is given by the subsidiary condition (j5.11|) .] It is straightforward to see from the Dirac 
equation for the monopole ()5.15p and the construction ()5.12|) that introducing the 
auxiliary dual field (which is a functional of F^jy and depends on the string function 
/.) 

1 

47r 

results in the following Dirac equation for the monopole field 



B,{x) = -— / {dy)r {x - y) *F,,{y) , (5.16) 



{i'-yd + gjBi^x) — m^) x{x) = 0. (5-17) 

Here we have chosen the string to satisfy the oddness condition [this is the "symmetric" 
solution, generalizing ()3.23|) ] 

rix) = -n-x), (5.18) 



which as we have seen is related to Schwinger's integer quantization condition 
Now ()5.14|) and ()5.17p display the dual symmetry expressed in Maxwell's equations ()5.2|] 
and ()5.4j) . Noting that satisfies [like taking Am = in (j3.16&j) 

(dx')r(x - x')B^{x') = 0, (5.19) 

we see that ()5.1(j|l is a gauge-fixed vector field [HZl IHHl defined in terms of the field 
strength through an inversion formula (see section 15.4.11) . In terms of these fields the 
"dual-potential" action can be re-expressed in terms of the vector potential and field 
strength tensor F^,^ (where B^ is the functional ()5.16|) of -F^,y) in first-order formalism 

as 

W = [ (dx)| - ^F^'^(x) {d^A, (x) - d,A^ {x)) + -^F^,{x)F^^{x) 

-|- ip^x) {i'-yd + 67^4(0;) — m^) ip^x) + x{x) {i'yd + g-yB{x) — m^) x{x)^, (5.20a) 
or in terms of dual variables, 

W = J (dx)| - ^ *F^''ix) (d^B, (x) - d,B, (x)) + ^ ^F^^ix) *F^,{x) 

+ il>{x) (i7(9 + 67^4(0;) — m^) il){x) + x{x) {}ld + 9lB{x) — m^) x(s) L (5.20&) 
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In (j5.2(Jaj) . A^{x) and F^^{x) are the independent field variables and B^{x) is given by 
()5.16p . while in ()5.20?)|) the dual fields are the independent variables, in which case, 

A,{x) = -^J idy)r {x - y) = ^^,uxa J {dy)nx - y) *F'^{y). (5.21) 

[Note that (jOM may be obtained from the form ()5.20aj) by inserting ()5.2H) into the 
former and then identifying according to the construction ()5.16|) . In this way the 
sign of -j^F^^Ff^^ = —^*F^^*F'^^ is fiipped.] Consequently, the field equation relating 
*F^''' and B^' is 

*F^, = d,B, - d,B, - I (dy) * {f,{x - y)jM - Ux - y)j,{y)) , (5.22) 
which is simply obtained ()5.10p by making the duality transformation ()2.2a|) . 

5.4- Quantization of dual QED: Schwinger-Dyson equations 

Although the various actions describing the interactions of point electric and magnetic 
poles can be described in terms of a set of Feynman rules which one conventionally uses 
in perturbative calculations, the large value of ag or eg renders them useless for this 
purpose. In addition, calculations of physical processes using the perturbative approach 
from string-dependent actions such as ()5.2()aj) and ()5.20^|1 have led only to string 
dependent results jHH]- In conjunction with a nonperturbative functional approach, 
however, the Feynman rules serve to elucidate the electron-monopole interactions. We 
express these interactions in terms of the "dual-potential" formalism as a quantum 
generalization of the relativistic classical theory of section 15.31 We use the Schwinger 
action principle [HIH 1^ to quantize the electron-monopole system by solving the 
corresponding Schwinger-Dyson equations for the generating functional. Using a 
functional Fourier transform of this generating functional in terms of a path integral 
for the electron-monopole system, we rearrange the generating functional into a form 
that is well-suited for the purpose of nonperturbative calculations. 



5.4-1 ■ Gauge symmetry In order to construct the generating functional for Green's 
functions in the electron-monopole system we must restrict the gauge freedom resulting 
from the local gauge invariance of the action ()5.20a|l . The inversion formulae for 
and i?^, ()5.21|1 and ()5.1(i|l respectively, might suggest using the technique of gauge-fixed 
fields IHZl as was adopted in jSH]. However, we use the technique of gauge fixing 
according to methods outlined by Zumino [HSl IMj and generalized by Zinn- Justin [HS] 
in the language of stochastic quantization. 

The gauge fields are obtained in terms of the string and the gauge invariant field 
strength, by contracting the field strength ()5.10|) . ()5.12|) with the Dirac string, f^{x), 
in conjunction with ()5.7|) . yielding the following inversion formula for the equation of 
motion, 

A^{x) = [ {dx')r{x - x')F^,{x') + d^K{x), (5.23) 
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where we use the suggestive notation, Ae{x) 

Ae(x) = ^ [ (dx')r(x - x')A, (x') . (5.24) 



4:71 _ 

In a similar manner, given the dual field strength (j5.22|) the dual vector potential takes 
the following form [cf. ^J^, (ICT^ ] 

B^{x) = -^J {dx')r{x - x') *F^,{x') + d^Ag, (5.25a) 

where 

Agi^) = {dx')nx - x')B,{x'). (5.256) 

It is evident that ()5.23|) transforms consistently under a gauge transformation 

A,{x) A,{x) + d,A,{x), (5.26) 

while in addition we note that the Lagrangian (|5.2(Jap is invariant under the gauge 
transformation, 

ip exp [ieAe] ^, A^-^ Af, + d^A^, (5.27a) 

as is the dual action (j5.20?)|l under 

X exp [i^^A^] X, B^^B^ + d^Ag. (5.276) 

Assuming the freedom to choose Ae(x) = — Ae(x), we bring the vector potential into 
gauge-fixed form, coinciding with 1)5.211) . 

where the gauge choice is equivalent to a string- gauge condition 

(da;')r(a; - x')A^{x') = 0. (5.29) 



[This is the analog of ()5.19p . and is equivalent to the gauge choice Ag = 0, see ()3.16aj) . 
used in section l!^.1.3l It is worth noting the similarity of this condition to the Schwinger- 
Fock gauge in ordinary QED, x ■ A{x) = which yields the gauge-fixed photon field, 
Afj,{x) = —x''jQdssF^iy{xs).] Taking the divergence of ()5.28j) and using ()5.2|) . the 
gauge-fixed condition ()5.28|) can be written as 



d,A'^ = J idy)rix-y)Uy), (5.30) 

which is nothing other than the gauge-fixed condition of Zwanziger in the two-potential 
formalism |8Uj . 

More generally, the fact that a gauge function exists, such that Ae(x) = — Ae(x) 
[cf. ()5.24|) ]. implying that we have the freedom to consistently fix the gauge, is in fact 
not a trivial claim. If this were not true, it would certainly derail the consistency of 
incorporating monopoles into QED while utilizing the Dirac string formalism. On the 
contrary, the string-gauge condition ()5.29|1 is in fact a class of possible consistent gauge 
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conditions characterized by the symbohc operator function (j5.8p depending on a unit 
vector ra'^ (which may be either spacehke or timehke). 

In order to quantize this system we must divide out the equivalence class of field 
values defined by a gauge trajectory in field space; in this sense the gauge condition 
restricts the vector potential to a hypersurface of field space which is embodied in the 
generalization of ()5.29|) 



^ I idx')nx~x')A,{x')=Ai{x), 



(5.31) 



where here A{ is any function defining a unique gauge fixing hypersurface in field space. 

In a path integral formalism, we enforce the condition (j5.31|) by introducing a 5 
function, symbolically written as 



5 (^/"^/^ - A{) = J [dAe] exp 



(dx)Ae(x) 



x( — I {dx')nx-x')A^{x')-Ai{x) 



(5.32) 



or by introducing a Gaussian functional integral 



[dXe] exp 



- I {dx){dx')\e{x)M{x,x')\e{x') 



+ i 



1 



(dx)Ae(a;) ( — / {dx')nx - x')A^{x') - A{(x) 



(5.33) 



where the symmetric matrix M{x, x') = k,^^6{x — x') describes the spread of the integral 
/ {dx')f^{x — x')A^{x') about the gauge function, A{(x). That is, we enforce the gauge 
fixing condition ()5.3H) by adding the quadratic form appearing here to the action ()5.20aj] 
and in turn eliminating Ae by its "equation of motion" 

(5.34) 



Ae(x) = ^[-^j (dy) Fix - y)A,{y) - A{(x) ) . 
Now the equations of motion ()5.2|) take the form, 

{dx')Xe{x')f^{x' -x) = 47rj^(a;) 



d'F^^ix) 



*F^,{x) 



{dx')\g{x')f^{x' -x) = 47r%(x) 



(5.35a) 
(5.355) 



where the second equation refers to a similar gauge fixing in the dual sector. Taking the 
divergence of (j5.35aj) implies Ae = from (j5.7|) and (j5.4j) . which consistently yields the 
gauge condition ()5.3H) . Using our freedom to make a transformation to the gauge-fixed 
condition ()5.28|) . A{ = 0, the equation of motion ()5.35aj) for the potential becomes 

A'^ix) = Anj^ix) + e.^^rjz-^d'^Yix), (5.36) 



{n ■ d) 



where we now have used the symbolic form of the string function ()5.8p . Even though 
fl5.34|) now implies n^^A^ = 0, we have retained the term proportional to n^n^, in the 
kernel, scaled by the arbitrary parameter k, 



K, 



(5.37) 
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so that i^ui, possesses an inverse 



9,,u —f + n 1 



D+{x), (5.38) 



n ■ d) \ Ann 11? J ■ Oy 

that is, J {dx')K^a{x — x')D°"^{x' — x") = g'^^5{x — x") , where -D+(x) is the massless scalar 
propagator, 

D+{x) = J . 5ix). (5.39) 

This in turn enables us to rewrite ()5.36p as an integral equation, expressing the vector 
potential in terms of the electron and monopole currents. 



A^{x) = 47r y {dx')D^,{x - x')f{x') 

+ e'^^- I {dx'){dx")D^,{x - x')h{x' - x")d';*Mx"). (5.40) 
The steps for Bfj,{x) are analogous. 

5.4-2. Vacuum persistence amplitude and the path integral Given the gauge-fixed 
but string-dependent action, we are prepared to quantize this theory of dual QED. 
Quantization using a path integral formulation of such a string-dependent action is by 
no means straightforward; therefore we will develop the generating functional making 
use of a functional approach. Using the quantum action principle (cf. 1,901121]) we write 
the generating functional for Green's functions (or the vacuum persistence amplitude) 
in the presence of external sources /C, 

Z{iq = (0+ I 0_)^, (5.41) 

for the electron-monopole system. Schwinger's action principle states that under an 
arbitrary variation 

5(0+ I 0_)'^ = i(0+ \5W{1C)\Q^)^ , (5.42) 

where W{'IC) is the action given in ()5.20aj) externally driven by the sources, /C, which 
for the present case are given by the set {J, *J, ?7, ^, 



W{1C) = W + j (dx) { JM^ + *J^^B^ + + + + , 



(5.43) 



V {v)j C iO being the sources for electrons (positrons) and monopoles (antimonopoles), 
respectively. The one-point functions are then given by 



i6J^'{x) ^ ' ' (0+|0_)'c ' i5*J^(x) ^ ' ' (0+|0_)^ ' 

S , 7.^. (0+1^(^)10-)"^ s (0+lx(^)|0-)^ ..... 

Using ()5.44j) we can write down derivatives with respect to the charges (here we redefine 
the electric and magnetic currents j ej and *j g*j) in terms of functional 
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5AJx) SJ>^{x) 



(0+I0-) 



-^{0+\0-)^ = i(0+| / idx)Tix)B^ix)\0^)'' = -i / (dx) 



6 



6 



6B„(x)6*J>^(x) 



(0+|0_)^. 
(5.45) 



Here we have introduced an effective source to bring down the electron and monopole 
currents, 



6 



SA,, 



i Srj 5fj' 



i S^^ <5f 



(5.46) 



•-11 "I "I ^--^11 
These first order differential equations can be integrated with the result 

6 6 



K. 



(o+|o_) 



— ie J (dx) 



exp 



ig j (dx) 



5A^{x) SJ>^{x) 



6B,{x)S*J'^{x) 
(0+|0-)^, 



(5.47) 



where (0+ | 0_)^ is the vacuum amplitude in the absence of interactions. By 
construction, the vacuum amplitude and Green's functions for the coupled problem 
are determined by functional derivatives with respect to the external sources /C of the 
uncoupled vacuum amplitude, where (0+ | 0_)q is the product of the separate amplitudes 
for the quantized electromagnetic and Dirac fields since they constitute completely 
independent systems in the absence of coupling, that is. 



(0+|0_)^ = (0+|A^(x)|0_)^ 



/c 



(5.48) 
(5.49) 



First we consider (0+ | 0_)q as a function of J and *J 
5 

i(5J'^(x) 

Taking the matrix element of the integral equation ()5.40|) but now with external sources 
rather than dynamical currents we find 



(0+|A^(x)|0_)^ = J {dx')D,,{x - x') \^A7rr{x') + e^'^^ J {dx")Mx' - x")cfi*Ux") 

x(0+|0_)^. (5.50) 
Using (j5.36p we arrive at the equivalent gauge-fixed functional equation. 



6 



■(0+|0-)^ 



47rt/^(x) -|- a^yuT 



i5J''{x) 

{dx')nx-x')d"^*r{x')] (0+|0_)^, 



which is subject to the gauge condition 



5J' 



(5.51) 



(5.52a) 
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or 



In turn, from ()5.47|) we obtain the full functional equation for (0+ | 0_)'^: 

6 



(5.526) 



9iiud + d^d^ + A'KKn^{n ■ d) 



^9 / (dy) 



exp 



X AttJJx) + e 



6 



6 



i6J''{x) 
- ie / (dy) 



(0+|0_)^ 
5 



fiuar 



idx')r{x - x')d"'*r{x') ) (0+ I 0_)^. 



(5.53) 



Commuting the external currents to the left of the exponential on the right side of ()5.53p 
and using ()5.44|) . we are led to the Schwinger-Dyson equation for the vacuum amplitude, 
where we have restored the meaning of the functional derivatives with respect to A, B 
given in ()5.46p . 

1 S 



i(5?7(x) i5r](x) 



-flUCTT 



iSJu{x) 



AttJ^{x) + 



idx')r{x - x')d"'*r{x') (0+ I 0_) 



(0+|0_)^ = (0+|i?,,(^)|0-)^, 



K 



In an analogous manner, using 
5 

i5*Ji'{x) 

we obtain the functional equation (which is consistent with duality) 

5 



(5.54) 



(5.55) 



i6*Ju(x) 



mx) "mx 



(dx')r(x -a:')5-e-|-7--47Y 1^(0+ I 0-) 
iori[x') 10 7] [X ) ' 



(dx')r(x - x')d"'.rix') (0+ I 0_) 



which is subject to the gauge condition 



{dx')f^'{x-x')- 



6 



In a straightforward manner we obtain the functional Dirac equations 



(5.56) 



(5.57) 



17^ + (77^- 



6 



i6J^'{x) 
6 

i6*Jt'(x) 



m 



6 



(0+ I 0_)^ = -r7(x)(0+ I 0_)^, (5.58a) 



i6ri{x) 
6 

i6^{x) 



(0+ I 0_)'^ = -^(a;)(0+ I 0_)^. (5.586) 



K 



In order to obtain a generating functional for Green's functions we must solve 
the set of equations (IH3T)|1 . (I5.58a|) . ()5.586jl subject to ()5.526jl and (pnTTfl for 
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(0+ I 0_)"^. In the absence of interactions, we can immediately integrate the Schwinger- 
Dyson equations; in particular, ()5.56|) then integrates to 

(0+1 O^)^'*-^ =Ar(J)exp|27ri J {dx){dx')%{x) D^'^ix - x')*J^ix') 

+ j (dx)(dx')(dx") V^(x)D^^(x - x')d'T{x' - x")J"(x")}. (5.59) 
We determine A/", which depends only on J, by inserting ()5.59p into ()5.54j) or ()5.51|) : 

InAr(J) = 27ri J {dx){dx')J^{x)D>"'{x - x')J^.{x'), (5.60) 
resulting in the generating functional for the photonic sector 
(0+ I 0-)lf''-^^ = exp {27ri [ {dx){dx')J^{x)D>"'{x - x')J^(x') 



+ 27ri J {dx){dx')*Jf,{x)D^''{x - x')*J^{x') 

-47ri [ {dx){dx')J^{x)D'"'{x - x')*J^{x")], (5.61) 



where we use the shorthand notation for the "dual propagator" that couples magnetic 
to electric charge 

D^. {x - x') = ^e^,,, j {dx")D+ {x - x") d"^ T {x" - x') . (5.62) 

The term coupling electric and magnetic sources has the same form as in ()5.5|) : here, 
we have replaced g'^'^D^, because of the appearance of the Levi-Civita symbol 

in ()5.62j) . [Of course, we may replace D^^'^ — > g^^D^ throughout ()5.6H) . because the 
external sources are conserved, d^J^ = d^*J^ = 0.] In an even more straightforward 
manner ()5.58a|l . (j5.58&jl integrate to 

(0+ I O^)^^'"'^'^^ = exp |i [ (dx)(dx') [f]{x)G4x - x')7]{x') + ^{x)G^{x - x')^{x')] 



(5.63) 

where and G^ are the free propagators for the electrically and magnetically charged 
fermions, respectively, 

GJx) = — 6{x), GJx) = — 6ix). (5.64) 

In the presence of interactions the coupled equations ()5.54|1 . ()5.56|1 . ()5.58ajl . (j5.58&j) are 
solved by substituting ()5.61|) and ()5.63p into ()5.47|) . The resulting generating function 
is 

Z(/C) = exp ( -ie / (dx)-^7'' 



6ri{x) i6J^^{x) 6ri{x) 
X exp I I (d,)^,'^^) Z„iK). (5.65) 
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5.4-3. Nonperturbative generating functional Due to the fact that any expansion in ag 
or eg is not practically useful we recast the generating functional ()5.65|) into a functional 
form better suited for a nonperturbative calculation of the four-point Green's function. 

First we utilize the well-known Gaussian combinatoric relation [Q^ llUUj : moving the 
exponentials containing the interaction vertices in terms of functional derivatives with 
respect to fermion sources past the free fermion propagators, we obtain (coordinate 
labels are now suppressed) 



Z{]C) = exp <^ i f]{G 



7] + Tr In ( 1 - 67 ■ r^G^ 



X exp < i / M G 



1-^7 



^ + Tr In 1-5-7 



i6*J 



Now, we re-express (j5.61|) . the noninteracting part of the generating functional of the 
photonic action, Zq{J, *J), using a functional Fourier transform. 



ZoU V) = / [dA] [dB] Zo (A, B) exp 



i / {J -A+V-B) 



or 



[dA] [d5]exp (iro[A,5,J, *J]) 



where (using a matrix notation for integration over coordinates) 



TT 



( J ■ A + V • 5) - — / A^K^.A'' + / B'^A-lB"" 



with the abbreviation 



1 



and the string-dependent "correlator" 
1 



{dx')dT{x - x')A^{x') 



(5.67a) 
(5.676) 

(5.68) 
(5.69) 



(47r)2 

Using ()5.68p we recast ()5.66|) as 



{dx") {rix - x")^ {x" - x') g,, -f,{x- x") U {x" - x')] .(5.70) 



Z(/C) 



[dA] [dB] FM)F2{B) exp (iro[A, 5, J, *J]) . 



(5.71) 



Here the fermion functionals Fi and F2 are obtained by the replacements ^ ^ A, 

5 



(5.72 a) 
(5.72&) 



F,{A) = exp jxr In (1 - 67 ■ AG^) + i j r] {G^[l - ej ■ AG^]'') r/| , 

F^iB) = exp jxr In {I - gj ■ BG^) {G^ [1 - 57 ■ BG^]'') . 

We perform a change of variables by shifting about the stationary configuration of the 
effective action, Vq[A,B, J, *J]: 

A^{x) = A^ix) + 0^(x), B'^ix) = B'^ix) + 0;(x) (5.73) 
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where A and B are given by the solutions to 

SFq {A,B,J,*J) ^ STq {A,B,J,*J) ^ 

6 A- ' 6B- ' ^ ^ ^ 

namely (most easily seen by regarding A and B' as independent variables), 

A^{x) = j {dx')D^,^{x - x') [A-Kr{x') - e"''"" j {dx")d'J^{x' - x") *J,(x")) , (5.75a) 

B^{x) = j {dx')D^,,{x - x')(^4nV^{x') + e^''''^ j {dx'')dXi^' - ^IM^'')), (5-756) 

reflecting the form of ()5.40|) and its dual. Note that the solutions ()5.75a|) . ()5.756|1 
respect the dual symmetry, which is not however manifest in the form of the effective 
action ()5.68|) . Using the properties of Volterra expansions for functionals and performing 
the resulting quadratic integration over (f){x) and (p'lx) we obtain a rearrangement 
of the generating functional for the monopole-electron system that is well suited for 
nonperturbative calculations: 

~ LJ iijj (X X ) 



exp <|27ri j (dx)(dx' 



Z^{l*J) "I J ' \6A^{x) '6A,{x') 

+ ^= , , D^^{x - x') - . . ) - 47ri / {dx){dx') - . . D^^ {x - x'^ 



5B^{x) "'"^ 6B^{x')J 6A^{x) ^ ' 6B,Ax' 



X exp |i J (d.)(d.O, (.) G(., + i / (d.)(d.')^ (.) .'m 

xexpj-^ de Ti-fAG{x,x\A) - dg' Ti -fBG{x, x\B)y (5.76) 

Here the two-point fermion Green's functions G{xi,yi\A), and G{x2,y2\B) in the 
background of the stationary photon field A, B are given by 

G{x,x'\A) = (x|(7p + -e-fAy^\x') , (5.77a) 
G{x,x'\B) = {x\{-fp + m^~ g-fBy^\x') , (5.776) 

where the trace includes integration over spacetime. This result is equivalent to the 
functional Fourier transform given in ()5.67a|) including the fermionic monopole-electron 
system: 

Z(/C) = j [dA] [dB] det (-17^'^ + m^) det (-17^5 + m^) 

xexpji J {dx){dx')(^f] {x)G{x,x'\A)t]{x') +^{x)G{x,x'\B)^{x')^^ 

X exp {-^ / {a^'K.^A'' - B'^K-lB"") + i j (J ■ A + *J ■ B)^ (5.78) 

where we have integrated over the fermion degrees of freedom. 

Finally, from our knowledge of the manner in which electric and magnetic 
charge couple to photons through Maxwell's equations we can immediately write the 
generalization of ()5.7(i|) for dyons, the different species of which are labeled by the index 
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a: 



Z{IC) = exp <j 2m / {dx){dx')J^'{x)V^^ {x - x') J" {x') 

-V.j Jx - X 



SA^ix) 'bA,{x' 



xexp|27ri ^ (da;)(da;'^ 

xexpji^ / (dx)(dx')Ca(a;)G'a(x,x'|A)Ca(a;') 

a 

xexpj-V/" dg Tr7AG'a(a;,x|gA)|. (5.79) 

where Aa = eaA + gaB, (a is the source for the dyon of species a, and a matrix notation 
is adopted, 

J^(x) - ( ^ - ( ^ (5 80a) 

and 

v., {x -x')=( ^""^ ["^ ~ ^2 - ) . (5.81a) 

' \ D^, {x - x') D^, {x-x') j ^ ' 

5.4-4- High energy scattering cross section In this subsection we provide evidence 
for the string independence of the dyon-dyon and charge-monopole (the latter being 
a special case of the former) scattering cross section. We will use the generating 
functional ()5.79|1 developed in the last subsection to calculate the scattering cross section 
nonperturbatively. We are not able in general to demonstrate the phenomenological 
string invariance of the scattering cross section. However, it appears that in much 
the same manner as the Coulomb phase arises as a soft effect in high energy charge 
scattering, the string dependence arises from the exchange of soft photons, and so 
in an appropriate eikonal approximation, the string-dependence appears only as an 
unobservable phase. 

To calculate the dyon-dyon scattering cross section we obtain the four-point Green's 
function for this process from ()5.79|) 

S 6 6 s 
G{xi,yi;x2,y2) = .^r , . Z{1C) . (5.82) 

i5Ci{xi) i5Ci(l/i) 1^(2(2:2) 1^(2(1/2) ^=0 

The subscripts on the sources refer to the two different dyons. 

Here we confront our calculational limits; these are not too dissimilar from those 

encountered in diffractive scattering or in the strong-coupling regime of QCD |1UH 

llU2t \10'6\ llU4t llU5j . As a first step in analyzing the string dependence of the 

scattering amplitudes, we study high-energy forward scattering processes where soft 

photon contributions dominate. In diagrammatic language, in this kinematic regime 

it is customary to restrict attention to that subclass in which there are no closed 

fermion loops and the photons are exchanged between fermions |lUlj . In the context 

of Schwinger- Dyson equations this amounts to quenched or ladder approximation (see 
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Figure 11. Dyon-dyon scattering amplitudes in the quenched approximation. 



figure ITT| . In this approximation the hnkage operators, L, connect two fermion 
propagators via photon exchange, as we read off from (j5.79j) : 

6 



exp{4../(d.)(d.0^2^-> 



(5.83) 



In this approximation ()5.82j) takes the form 



G {xi,yi;x2,y2) = -e "G'i(xi, A)G'2(a;2, Z/2I A 



(5.84) 

A=B=0 

where we express the two-point function using the proper-time parameter representation 
of an ordered exponential 



Ga{x,y\Aa) = i / d^e 



-i^ima-iyd) 



expji^ de'e«'^Mae-^'^^| 



5{x-y), 



(5.85) 



where denotes path ordering in The 12 subscripts in L12 emphasize that only 
photon lines that link the two fermion lines are being considered. 

Adapting techniques outlined in |106[llU7j we consider the connected form of ()5.84p . 
We use the connected two-point function and the identities 



and 



5Af,{x 



dae^^L 



-G{y,z\A) = eG{y,x\A)-f^G{x,z\A). 



(5.86) 



(5.87) 



Using ()5.84|) and ()5.85|) one straightforwardly is led to the following representation of 
the four-point Green function, 

G'(xi,?/i;X2,?/2) = -4vri / da {dzi){dz2) [qi ■ q2 D^^{zi - Z2) - qiXq2Df,^{zi - Z2) 



xe 



aLi 



Giixi, zi\Ai)j^Gi{zi, yi\Ai)G2{x2, z2\A2)YG2{z2, y2\A2] 



A=B=0 

where the charge combinations invariant under duality transformations are 

-m'hc = —Kc. 



qi ■ q2 = 6162 + 9192 = q, 



qixq2 = ei92 - 9162 



(5. 



(5.89) 
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In order to account for the soft nonperturbative effects of the interaction between 
electric and magnetic charges we consider the hmit in which the momentum exchanged 
by the photons is small compared to the mass of the fermions. This affords a substantial 
simplification in evaluating the path-ordered exponential in ()5.85|) : in conjunction 
with the assumption of small momentum transfer compared to the incident and 
outgoing momenta, q/p{i^2) ^ 1, this amounts to the Bloch-Nordsieck |lU8j or eikonal 
approximation (see [109, llUl 11111 11121 IllJl 11141 1115j : for more modern applications in 
diffractive and strong coupling QCD processes see |1()H ll()2| ll()3| ll()4| ll()5j ). In this 
approximation ()5.85p becomes 



(5.90) 



G,{x,y\A)^i rdee-«™5(x-i/-e-)exp(i / d^'^ ■ ^ (x - e'- 

With this simplification each propagator in ()5.84|) can be written as an exponential 
of a linear function of the gauge field. Performing mass shell amputation on each 
external coordinate and taking the Fourier transform of ()5.88|) we obtain the scattering 
amplitude, T{pi,p[;p2,P2)- 

T{pi,p^,^P2,p'2) ^ l^ae'^^^^l {dz,){dz2) (qi ■ D^,{z, - Z2) - q,xq2D,,{z, - Z2)) 
X / {dx,)e-'P^'-'^u{p,) (mi + V, ■ p,) G,{x,, z,\A,)r 



X j {dy^)e'^^y^G,{z,,yM^) {m, + v[ ■ p[) u{p[) 

X j {dX2)e-'^'^'u{p2) {m2 + V2 ■ P2) G2{X2, Z2\A2K 

X j {dy2y^'-y'G2{z2, y2\A2) (ms + v'^ ■ p',) u{p',). (5.91) 
Substituting ()5.90|) into ()5.91|) . we simplify this to 

T(pi,K;P2,P2) ^ /(d2i)(d22)e~'"^(^^~^'^)e-'"^(p^-^2)u(p;)7'^M(pi)M(p^)7"M(p2) 



-47ri 

X 



(qi ■ q2 D^^uizl - Z2) - qiXq2D^^{zi - 2:2)) e''^^^ 

POO 

i / dai {pi ■ Ai {zi + aipi) + p[ ■ Ai {zi - aip[) } 

POO 

i / da2 {p2 ■ A2 {z2 + 02^2) +P2- -^2 {Z2 - "2^2) } 

, ^0 



X exp 



X exp 



(5.92) 



Choosing the incoming momenta to be in the z direction, in the center of momentum 
frame, p^ = {Ei, 0, 0,p), P2 = {E2, 0, 0, —p), invoking the approximation of small recoil 
and passing the linkage operator through the exponentials containing the photon field, 
we find from ()5.92j) 

T{pi,p'i;P2,P2) „ 



-Am 



da / (d^i)(d^2)e-'^^(^^-^'Oe-'^^(^'^-pO^(p;)7^M(pi)t^(P2)7.t^(P2) 



X (qi ■ q2D^'' {zi - Z2) - qixq2^^" (zi - ^2)) e^«*(Pi'f2;.i-.2)^ (5 93) 
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where the "eikonal phase" integral is 

= 47rpiP2 / dai da2 Ui ■ q2 DkX - qixqa^KAj (^1 - ^2 + am - a2P2) ■ (5.94) 

We transform to the center of momentum coordinates, by decomposing the relative 
coordinate accordingly, 

{zi - Z2Y = - Tip'^ + rsp^, (5.95) 

where the Jacobian of the transformation is 

J = p^s (5.96) 

and s = —{pi+P2)^ is the square of the center of mass energy. Here we use the symmetric 
infinite string function, as discussed in section El which has the momentum-space form, 

Inserting the momentum-space representation of the propagator, and recalling ()5.f)2|l . 
we cast ()5.94|) into the form 



^{pi,P2;x) ^ Anp'^p^ / daida2 



X 



qi ■ q2 9KX - qixq2e«A<7r/i; ^ ; + 



(5.9^ 



2 \n ■ k — ie n ■ k + ie^ 

where we have introduced the standard infrared photon-mass regulator, /x^. The delta 
functions that result from performing the integrations over the parameters ai and ^2 in 
()5.98|) in the eikonal phase suggests the momentum decomposition 

k^ = k^ + XiCi + X2G2, where Xi = P2 ■ k, and X2 = Pi ■ k, (5.99) 

and the four- vector basis is given by 

e^ = ^fl,0,0,^) and = A, 0, 0, -^V (5.100) 

which have the following properties, in terms of the masses mi and m2 of the two dyons, 

1 m? 1 1 Pi ■ P2 / 
ei ■ ei = 5-, 62 ■ 62 = ^, and ei ■ 62 = (5.101) 

g pz g pz s p"^ 

The corresponding measure is 

(dk) = J-M2k^dAidA2, (5.102) 

in terms of the Jacobian in ()5.9(ij) . Using the definition of the M0ller amplitude, M{s, t), 
given by removing the momentum- conserving delta function, 

nPi,p[;P2,P2) = {2nY6^^\P-P')M{s,t), (5.103) 
we put fl5.93j) into the form 



M{s,t)^-iJ^ da J d'xxe-'^^-"^u(p;)7^M(pi)M(p'2)7'^M(p2)Ve'"*^^''^''"\ (5.104) 
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d^k^ dAi dAa 



27r5(Ai)27r(5(A2) 



X 



1 1 



n 



+ 



2 \n ■ k — ie n ■ k + ie 



(5.105) 



Here P = Pi + P2 and P' = p[ + P2, and q = Pi — p'l is the momentum transfer. The 
factor 

" 1 



exp(iripi ■ q - iT2P2 ■ q) = exp 



+^2) 



(5.106) 



has been omitted because it is unity in the eikonal hmit, and correspondingly, we have 
carried out the integrals on ti and T2. The eikonal phase (j5.98j) now takes the very 
similar form 



^{pi,P2;x) 



P1P2 

Py/S 



(5.107) 



Choosing a spacelike string in order to have a local interaction in momentum space, 
= (0, n), integrating over the coordinates Ai,A2, and introducing "proper-time" 
parameter representations of the propagators, we reduce ()5.107|) to 



^{pi,P2;x) 



d^k 



p^sj (27r)2 
n'^ d 



x< qi • qapi -^2 - ciiy«^2PiP2^fiuaT 



2i dn^ 



it(n-k+je) 



2qi ■ q2 ^ i^o (/U |x|) - qi^o^2(^zjkn^ — 



P\fs 



dl 
it 

^iro(/w|(x + tn)|) 



dl 



^it(n-k— je) 



It 



(5.10^ 



in terms of modified Bessel functions, where we have dropped the subscript _L. 

We perform the parameter integral over t in the limit of small photon mass /x^: 



z ■ (n X x) 

2 ^ ^ 



1 



dte-'l*l 



(t + fi ■ x)2 + 



arctan 



n ■ X 



n ■ X 



z ■ n X X 



(5.109) 



so the phase is 

$(pi,P2;x) ^ 2 <! qi ■ q2ln (/i|x|) - qixq2 arctan 



n ■ X 



z ■ (^n X X 

In this limit we have used the asymptotic behavior of the modified Bessel function 



(5.110) 



In 



X 0, 



(5.111) 



where 7 = 0.577... is Euler's constant and we have defined /i = e'^/i/ 2. Similarly, 
fl5.104p becomes 



M(s,t) ^ -2i / da / d^:i^e-''^■''u{p\)^^'u{pl)u{p'2)Yu{p2) 



y-lgfiu^i ■q2^o(Ai|x|) - e^^^^qixq2ra' 



d 1 
9^2 



dl 



iro(/i|(x + tn) 



(5.112) 
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Although in the eikonal hmit, no spin-flip processes occur, it is, as always, easier 
to calculate the helicity amplitudes, of which there is only one in this case. In the 
high-energy limit, p° ^ m, the Dirac spinor in the helicity basis is 



w"(p) = '^|^(l + i75aK, (5.113) 

where the may be thought of as two-component spinors satisfying 7°Wo- = v^- They 
are further eigenstates of the helicity operator cr ■ p with eigenvalue a: 

vlip') = (^cos ^, sin , vl{p') = sin ^, cos , (5.114a) 

MP) = ( J ) ' ^-(P) = ( 1 ) • ^^-^^^^^ 

We employ the definition 

75 = 7°7SV (5.115) 

and consequently 7^7 = i7^cr, where aij = eijkcr^. We then easily find upon integrating 
over the parameter a that the spin nonflip part of ()5.112j) becomes {6 — > 0) 

M(s,t) = 1 [ d^^e-i^--emPi,P2;-) _ (27i?S\q)\. (5.116) 

2mim2 t J J 

Now notice that the arctangent function in ()5.110|) is discontinuous when the xy 
component of n and x lie in the same direction. We require that the eikonal phase factor 
e** be continuous, which leads to the Schwinger quantization condition (|3.2|) : 

qixq2 = -m', (5.117) 

where m' is an integer. Now using the integral form for the Bessel function of order 1/ 



2n 



27r 







we find the dyon-dyon scattering amplitude ()5.11fj|l to be [see also ()7.13|1 below] 

M(s,t) = — ^—e-^^rn'^ dxxJ2.„,(gx)e'2"i°(^"), (5.119) 
mima Jo 

where a = qi ■ q2, and ip is the angle between qj_ and fij,. The integral over x is just a 
ratio of gamma functions. 

Then ()5.119|1 becomes 

M(,, i) « -^^(,n' - ia)e--'* ( '4)'^ + + (5.121) 



17111712 \ J V {1 + m' — \a) 

This result is almost identical in structure to the nonrelativistic form of the scattering 
amplitude for the Coulomb potential, which result is recovered by setting m' = 0. (See, 
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for example, Following the standard convention |117j we calculate the spin- 

averaged cross section for dyon-dyon scattering in the high energy limit, 



While the Lagrangian is string-dependent, because of the charge quantization condition, 
the cross section, ()5.122j) . is string independent. Not surprisingly, this coincides with 
the Rutherford formula (jTTHIl . (j;-i.i;-i5jl . 

For the case of charge-monopole scattering ei = g2 = 0, this result, of course, 
coincides with that found by Urrutia jll8j . which is also string independent as a 
consequence of fl3.1|) . We should also mention the slightly earlier work of Ore, 
demonstrating the Lorentz invariance of charge-monopole scattering |119j . This is to be 
contrasted with ad hoc prescriptions that average over string directions or eliminate its 
dependence by simply dropping string-dependent terms because they cannot contribute 
to any gauge invariant quantities (cf. [5^). 

5.5. Conclusion 

In this section we have responded to the challenge of Schwinger [SS], to construct a 
realistic theory of relativistic magnetic charges. He sketched such a development in 
source theory language, but restricted his consideration to classical point particles, 
explicitly leaving the details to the reader. Urrutia applied this skeletal formulation 
in the eikonal limit |118j . as already suggested by Schwinger. 

We believe that we have given a complete formulation, in modern quantum field 
theoretic language, of an interacting electron-monopole or dyon-dyon system. The 
resulting Schwinger-Dyson equations, although to some extent implicit in the work of 
Schwinger and others, were given in [HI] for the first time. 

The challenge remaining is to apply these equations to the calculation of monopole 
and dyon processes. Perturbation theory is useless, not only because of the strength 
of the coupling, but more essentially because the graphs are fatally string- or gauge- 
dependent. The most obvious nonperturbative technique for transcending these 
limitations in scattering processes lies in the high energy regime where the eikonal 
approximation is applicable; in that limit, our formalism generalizes the lowest-order 
result of Urrutia and charts the way to include systematic corrections. More problematic 
is the treatment of monopole production processes — we must defer that discussion 
to subsequent publications. In addition we have also detailed how the Dirac string 
dependence disappears from physical quantities. It is by no means a result of string 
averaging or a result of dropping string-dependent terms as in [HH]- In fact, it is a 
result of summing the soft contributions to the dyon-dyon or charge-monopole process. 
There are good reasons to believe that inclusion of hard scattering contributions will 
not spoil this consistency. At the level of the eikonal approximation and its corrections 
one might suspect the occurrence of a factorization of hard string-independent and soft 
string-dependent contributions in a manner similar to that argued in strong-coupling 




(5.122) 
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QCD. 

It is also of interest to investigate other nonperturbative methods of calculation in 
order to demonstrate gauge covariance of Green's functions and scattering amplitudes 
in both electron-monopole and dyon-dyon scattering and in Drell-Yan production 
processes. In addition there is a formalism employed in |12()|ll2l]ll22j based on Fradkin's 
|123j Green's function representation, which includes approximate vertex and self-energy 
polarization corrections using nonperturbative techniques, which we are adapting to the 
magnetic charge domain. (For a first pedagogical example of this formalism see |124j .) 
We hope in the future to apply the techniques and results found here to the Drell-Yan 
production mechanism, for example, and obtain phenomenologically relevant estimates 
for the laboratory production of monopole-antimonopole pairs. 



6. Renormalization 

As discussed in sectional Lorentz invariance (rotational invariance for the nonrelativistic 
theory) is satisfied by the dual electrodynamics of electric and magnetic charges 
interacting provided the quantization condition is obeyed. But is the theory 
renormalizable? This question was addressed by Schwinger in 1966 [THllZZj. His view 
at the time was that renormalization described the connection between the particle and 
field level description of reality. At both of these levels consistency demanded that the 
quantization condition must hold, so that if integer quantization is appropriate, 

^— = ^0, ^ = ^' (6-1) 

he he 

but that the integers n and uq need not be the same. The question is whether electric 
and magnetic charges are renormalized by the same, or different factors. He argued that 
the former was the case, because charge renormalization refers to the electromagnetic 
field, not its sources. That is 

- = ^ = C < 1, (6.2) 

eo 9o 

so in view of the charge quantization condition ()6.1|) the quantum numbers no and n 
are not the same: 

= -. (6.3) 

no 

The discreteness of renormalization of the dual theory is thus manifest from this point 
of view. 

This is at odds with the modern understanding of renormalization as a continuous 
evolution of parameters, such as the charge, with change of energy scale. It would 
seem that this view of the renormalization group may be difficult to maintain without 
a perturbative framework: That is, at any energy scale Q, we might expect 

e{Q)g{Q) = n. (6.4) 

For this reason Laperashvili and Nielsen |125t I126t I127| I128| I129j , following Zwanziger 
[Hm IS21 argue that ()6.4j) holds at all scales, or in terms of the bare and renormalized 
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quantization numbers, n = Uq. That is, the electric and magnetic charges are 
renormahzed by exactly inverse factors. In terms of the fine structure constants, for 
the minimal Dirac pole strength, m' = 1/2, this says 

*a{Q)a{Q) = ^. (6.5) 

Laperashvili, Nielsen, and collaborators have exploited the small window which this 
seems to permit for perturbative calculations, where neither a nor *a are bigger than 
unity. 

However, at best there is room for serious doubt about the essential validity 
of this procedure. In ordinary quantum electrodynamics charge renormalization can 
be regarded as arising entirely from vacuum polarization. Presumably, this is still 
the case in dual QED. Using lowest order perturbative graphs to describe vacuum 
polarization does violence to the charge quantization condition; moreover, higher order 
graphs involving both electrically and magnetically charged particles necessarily bring 
in the Dirac string, which as we have repeatedly emphasized can only disappear in a 
nonperturbative treatment. Such is as yet lacking in our analysis of renormalization in 
dual QED. 

7. Eikonal approximation 

It is envisaged that if monopoles are sufficiently light, they would be produced by 
a Drell-Yan type of process occurring in pp collisions at the Tevatron. Two photon 
production channels may also be important. The difficulty is to make a believable 
estimate of the elementary process gg ^ 7* — MM, where q stands for quark and M for 
magnetic monopole. It is not known how to calculate such a process using perturbation 
theory; indeed, perturbation theory is inapplicable to monopole processes because of 
the quantization condition 1)3.11) . It is only because of that consistency condition that 
the Dirac string, for example, disappears from the result. 

Only formally has it been shown that the quantum field theory of electric and 
magnetic charges is independent of the string orientation, or, more generally, is gauge 
and Lorentz invariant [7^1 EHl |77| EHl E3 EHl EOl EH E21 • It has not yet proved possible to 
develop generally consistent schemes for calculating processes involving real or virtual 
magnetically charged particles. Partly this is because a sufficiently general field theoretic 
formulation has not yet been given; a small step in remedying this defect was given 
in jHlj, reviewed in sectional However, the nonrelativistic scattering of magnetically 
charged particles is well understood, as described in section El Thus it should not be 
surprising that an eikonal approximation gives a string-independent result for electron- 
monopole scattering provided the condition (jH.lj) is satisfied. In section El we described 
the eikonal approximation in terms of the full field-theoretic formulation. Since that 
formalism is rather elaborate, we give here a simplified pedagogical treatment, as 
described in |13Uj . and first worked out by Urrutia |118j . 
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The interaction between electric {J'^) and magnetic {*J^) currents is given by (j5.5p . 

or 

Wi"> ^ -._/(d.),d.')(d.'V(.)a'I>.(. - - (7.1) 

Here is the usual photon propagator, and the arbitrary "string" function /^(x — x') 
satisfies ()5.7j) . or 

djf'ix-x') =47r6{x-x'). (7.2) 

It turns out to be convenient for this calculation to choose a symmetrical string, which 
satisfies ()5.18|) . or 

Fix) = ~n-x). (7.3) 

In the following we choose a string lying along the straight line n'^, in which case the 
function may be written as a Fourier transform ()5.97|) . or 



2i J {2tt)^ \n ■ k — le n-k + iej 
In the high-energy, low-momentum-transfer regime, the scattering amplitude 
between electron and monopole is obtained from ()7.1|) by inserting the classical currents, 

r{x) = e r dX^6(x- ^X) , (7.5a) 
m \ m / 

*r{x) =g j^J\'^6[x + h- ^A') , (7.5&) 

where m and M are the masses of the electron and monopole, respectively. Let us 
choose a coordinate system such that the incident ultrarelativistic momenta of the two 
particles have spatial components along the 2;-axis: 

Pi = (p, 0,0,p), p2 = (p, 0,0,-p), (7.6a) 

and the impact parameter lies in the xy plane: 

6=(0,b,0). (7.66) 

Apart from kinematical factors, the scattering amplitude is simply the transverse Fourier 
transform of the eikonal phase, which is the content of ()5.116|) . 



/(q) = j d^be-^^ 'i (e'^-1). 



(7.7) 



where x is simply W^^^^ with the classical currents substituted, and q is the momentum 
transfer. 

First we calculate x'l it is immediately seen to be, if has no time component, 

d^ki z ■ fn X ki ) u / 1 1 



X = 27reg / J ^ ^^'^ [ ^ , ^ ■ + ^ .\ ■ ) ^ (7- 
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where is the component of the photon momentum perpendicular to the z axis. From 
this expression we see that the result is independent of the angle n makes with the z 
axis. We next use proper-time representations for the denominators in ()7.8|1 . 



oo 



^= / dse-^'^', (7.9a) 







1 1 

+ 



oo 

Jo J -oo 



, , e"^e-'^l^ (7.96) 

n ■ kx — ie n ■ k^ + ie i 

We then complete the square in the exponential and perform the Gaussian integration 
to obtain 

X = e^z ■ (n X b) / dA (7.10a) 

J_^ (A + b ■ n)^ + 6^ - (b ■ n)^ 



or 



Y = 2eg arctan I 3: — ^ rr I , (7.106) 

Vz ■ (b X n)y ^ ' 

which is contained in ()5.110|) . Because e'^ must be continuous when and b lie in 
the same direction, we must have the Schwinger quantization condition for an infinite 
string, 

eg = m' , (7-11) 

where m' is an integer. 

To carry out the integration in ()7.7|) . choose b to make an angle with q_|_, and 
the projection of n in the xy plane to make an angle with q_L; then 

X = 2eg{<j)-^-7,/2). (7.12) 

To avoid the appearance of a Bessel function, as occurs in ()5.119|) . we first integrate 
over b = |b|, and then over 0: 

/(q) = / d0 / l)^l)Q--M{cos<l,-ie)^2im'(^-^-^/2) 

Jo Jo 

4 -2im'(i/)+7r/2) f j 2m'-l ATrm' 

I az Z ^'^ "\^-2im'ii (7 13) 



i J(j {z + 1/ z — ieY (f 

where C is a unit circle about the origin, and where again the quantization condition 
()7.11|) has been used. Squaring this and putting in the kinematical factors we obtain 
Urrutia's result [mi [cf. 

% = 4vr(e,)^i, t = q\ (7.14) 

which is exactly the same as the nonrelativistic, small angle result found, for example, 
in (ITTH|1 . 
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8. Energy loss by magnetic monopoles traversing matter 

The essence of the energy loss mechanism of charged particles traveling through matter 
can be described by classical electrodynamics. By a simple duality analysis, therefore, 
one should be able to describe the rate at which a particle carrying magnetic charge 
loses energy when it passes through matter. The details of this argument can be found 
in the last two chapters of j2ZI- It is based on the fundamental analyticity requirements 
of the electrical permittivity, demanded by causality, the Kramers-Kronig relations. In 
terms of positive spectral functions p{uj) and ^(ti;), which satisfy 



du'piiu') 



the dielectric function obeys 



1 



1 + 



duj'- 



12 



doj'- 



-{uj + iey 



ijj 



12 



[oj + ie 



Here Up is the plasma frequency. 



OJr, 



m 



(8.1) 

(8.2a) 
{8.2b) 

(8.3) 



in terms of the electron mass m, and density of free electrons n. 

Using ()8.26|) . in Chapter 52 of |2Z] we derive the following formula for the energy 
loss —dE when a charged particle (charge Ze) having velocity v travels a distance dz: 



dE looliZef 



In- 



2„,2 



dz 2 t;2 

where the last integral should be omitted iiv/c< 1 
duj q{uj) Incj^ = IntUg, 



l/.(0) 



dr-Y^ 



e(0). Here 



•4) 



(8.5) 



and is given by the root of 1 — f^£:(iz/)/c^, that is, 

,2 



duj'- 



q[uj 



(8.6) 

Further, K is boundary between low momentum transfer events and high momentum 
ones, such as 5-rays. For a more complete theory, and extensive comparison with 
experiment, the reader is referred to |131j . Passage of particles through matter. 

In parallel with the above derivation, we can use ()8.2a|) to derive the corresponding 
formula for the energy loss rate by a magnetically charged particle: 

;^2^2 pe{0) / ^2 



dE 
dl 



2 c2 



In 



< (1 - 



1 



where the latter integral only appears if ^ > 

POD 

/ duj p{iu) \n Lu"^ = In uj^. 
Jo 



/e(0) 



=. Here 



(8.7) 



18) 
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We note that the predominant change from the energy loss for electrically charged 
particles lies in the replacement 
Ze q 

. ^, 8.9 

V c 

as earlier claimed in ()2.16|) . provided ujI ~ cu^. 

This result should only be regarded as quahtative. For our experimental analysis, 
we will use the extensive results of Ahlen |132t I133j . As we discussed in section 
Kazama, Yang, and Goldhaber jHSl have obtained the relativistic differential scattering 
cross section for an electron moving in the magnetic field of a fixed magnetic pole. Ahlen 
then used this cross section to obtain the following expression for monopole stopping 
power: 

where A^e is the number density of electrons, I is the mean ionization energy, = 
0.406 (0.346) is the Kazama, Yang and Goldhaber correction for magnetic charge 
2\m'\ = n = 1 {n > 2) respectively, S is the usual density correction and -B(|?t-|) = 0.248 
(0.672, 1.022, 1.685) is the Bloch correction for n = 1 {n = 2,3,6), respectively |131j . 
(Of course, one must divide by the density to get dx in g/cm^.) This formula is good 
only for velocities [3 = v/c > 0.1. For velocities [3 < 0.01, we use (60) of |132j as an 
approximation for all materials: 
(\E 

- = (45 Gev/cm)n2/5, (8.11) 
dx 

which is linear in j3 in this region. The two dE/dx velocity regions are joined by an 
empirically fitted polynomial in the region oi (3 = 0.01-0.1 in order to have a smooth 
function of (3. For the elemental and composite materials found in the DO and CDF 
detectors, we show the resulting dE/dx curves we used in figure IT^ (See |134j .) 

9. Binding 

To this point in this review, we have concentrated on the scattering of monopoles with 
charged particles, or on dyon-dyon scattering. Now we turn to the question of the 
binding of these particles. Our discussion will be largely, although not exclusively, 
based on the nonrelativistic description. 

If g = 6162 + gig2 < 0, and m' = —{eig2 — e2gi)/hc, TYnr (I3.126|) gives binding: 

(9.1) 



Enj = 



1 . . ,_.o ,2\l/2 



N + -+{{j + l/2)'-m") 



where A^ is a principal quantum number. We will not further address the issue of dyons 
|135| 11361 1137t I138t I139t I140j , which for the correct sign of the electric charge will always 
bind electrically to nuclei. Monopoles will not bind this way; rather, a magnetic moment 
coupling as in ()3.137j) is required; for example, for spin-1/2, 

^5 = -^7'^-B, 7 = 1 + «: = |. (9.2) 
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Figure 12. Energy loss of a magnetic monopole in various materials. These dE/dx 
curves are for a magnetic charge value of 2|m'| = n = 1; apart from the correction 
terms iir(|n|) and we multiply by for larger magnetic charge values. 



(7 = 1 or (yf = 2 is the "normal" value.) 

Suppose monopoles are produced in a collision at the Tevatron, for example; they 
travel through the detector, losing energy in a well-known manner (see, e.g., which 
results are summarized in section |H1), presumably ranging out, and eventually binding 
to matter in the detector (Be, Al, Pb, for example). The purpose of this section is to 
review the theory of the binding of magnetic charges to matter. 

We consider the binding of a monopole of magnetic charge g to a. nucleus of charge 
Ze, mass Ai = ArUp, and magnetic moment 

M = ^7S, (9.3) 



rUpC 



S being the spin of the nucleus. (We will assume here that the monopole mass ^ A4, 
which restriction could be easily removed.) The charge quantization condition is given 
by (jH.lj) . Because the nuclear charge is Ze, the relevant angular momentum quantum 
number is (m' being an integer or an integer plus 1/2) 

I = \m'\Z. (9.4) 
9.1. N onrelativistic binding for S = 1/2 

In this subsection we follow the early work of Malkus |141j and the more recent papers 
of Bracci and Fiorentini |142[ 11431 1144| H^H] . (There are also the results given in |14(ij . 
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but this reference seems to contain errors.) 

The neutron {Z = 0) is a special case. Binding will occur in the lowest angular 
momentum state, J = 1/2, if 

ItI > TT^ (9.5) 

Since 7„ = —1.91, this condition is satisfied for all m' . 

In general, it is convenient to define a reduced gyromagnetic ratio, 

7 = — 7, K = 7-l. (9.6) 

This expresses the magnetic moment in terms of the mass and charge of the nucleus. 
Binding will occur in the special lowest angular momentum state J = / — ^ if 

7 > 1 + ^. (9.7) 

Thus binding can occur here only if the anomalous magnetic moment k > 1/4/. The 
proton, with n = 1.79, will bind. 

Binding can occur in higher angular momentum states J if and only if 

\k\> n^ = ^\j^ + J (9.8) 

For example, for J = / + i, Kc = 2 + 3/4/, and for J = / + |, Kc = 4 + 15/4/. Thus 
^He, which is spin 1/2, will bind in the first excited angular momentum state because 
k = -4.2. 

Unfortunately, to calculate the binding energy, one must regulate the potential at 
r = 0. The results shown in tabled assume a hard core. 

9.2. Nonrelativistic binding for general S 

The reference here is jl47j . The assumption made here is that I > S. (There are only 
3 exceptions, apparently: ^H, *Li, and ^°B.) 

Binding in the lowest angular momentum state J = / — S* is given by the same 
criterion ()9.7|) as in spin 1/2. Binding in the next state, with J = I — S + 1, occurs if 
X± > J where 



A± = - |/ - 2/ - 1 ± y (1 + /)2 + (25 - 1 - /)|/ + \P (^l) (9.9) 

The previous result for S' = 1/2 is recovered, of course. = 1 is a special case: Then 
A_ is always negative, while A+ > ;^ if 7 > 7c, where 

^ 3 (3 + 16/ + 16/^) 
^^4/ 9 + 4/ ■ ^ 
For higher spins, both A± can exceed 1/4: 

A+ > i for 7 > 7c_ (9.11a) 
A_ > i for 7 > 7c+ (9.116) 
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where for S 



3 
2 



(9.12) 



For 4Be, for which 7 = —2.66, we cannot have binding because 3 > 7c- > 1.557, 
3 < 7c+ < 8.943, where the ranges come from considering different values of 2\m'\ from 
1 to 00. For 5* = |, 



So ^gAl will bind in either of these states, or the lowest angular momentum state, because 
7 = 7.56, and 1.67 > 7^- > 1.374, 1.67 < 7^+ < 4.216. 

9.3. Relativistic spin- 1/2 

Kazama and Yang treated the Dirac equation . See also |148[ 1149] and |135| I136| 

iTTfirniKirnMrran] . 

In addition to the bound states found nonrelativistically, deeply bound states, with 
-^binding = -M. are found. These states always exist for J > / + 1/2. For J = / — 1/2, 
these (relativistic) E = bound states exist only if k > 0. Thus (modulo the question 
of form factors) Kazama and Yang [HHl expect that electrons can bind to monopoles. 
(We suspect that one must take the existence of these deeply bound states with a fair 
degree of skepticism. See also |150j .) 

As expected, for J = / — 1/2 we have weakly bound states only for k > 1/4/, which 
is the same as the nonrelativistic condition ()9.7|) . and for J > / + 1/2, only if \k\ > k^, 
where Hc is given in ()9.8j) . 

9.4. Relativistic spin-1 

Olsen, Osland, and Wu considered this situation |15H I152j . 

In this case, no bound states exist, unless an additional interaction is introduced 
(this is similar to what happens nonrelativistically, because of the bad behavior of the 
Hamiltonian at the origin). Bound states are found if an "induced magnetization" 
interaction (quadratic in the magnetic field) is introduced. Binding is then found for 
the lowest angular momentum state J = / — 1 again if k > 1/4/. For the higher angular 
momentum states, the situation is more complicated: 

• for J = /: bound states require / > 16, and 

• for J > / + 1: bound states require J( J + 1) — /^ > 25. 

But these results are probably highly dependent on the form of the additional 
interaction. The binding energies found are inversely proportional to the strength A 
of this extra interaction. 



(7c) 



36 + 28/ T V1161 + 1296/ + 64/2 
12/ 



(9.13) 
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Table 1. Weakly bound states of nuclei to a magnetic monopole. The angular 
momentum quantum number J of the lowest bound state is indicated. In Notes, 
NR means nonrelativistic and R relativistic calculations; he indicates an additional 
hard core interaction is assumed, while FF signifies use of a form factor. IM represents 
induced magnetization, the additional interaction employed for the relativistic spin-1 
calculation. We use \m'\ = 1/2 except for the deuteron, where |rn'| = 1 is required for 
binding. 



Nucleus 


Spin 


7 


7 


J 








Notes 


Ref 


n 


1 

2 


-1.91 




1 

2 






350 keV 


NR,hc 


[nnj 




1 
2 


2.79 


2.79 


/ - 


1 

2 


= 


15.1 keV 
320 keV 
50-1000 keV 
263 keV 


NR,hc 
NR,hc 
NR,FF 
R 


[14S] 


?H 


1 


0.857 


1.71 


l- 


1 


= (|m'| = 1) 


if keV 


R,IM 


jTKTllTH^ 


iHe 


1 
2 


-2.13 


-3.20 


1 + 


1 
2 


_ 3 
2 


13.4 keV 


NR,hc 


m 


?^A1 


5 
2 


3.63 


7.56 


I - 


5 
2 


= 4 


2.6 MeV 


NR,FF 




?^A1 


5 
2 


3.63 


7.56 


/ - 


5 
2 


= 4 


560 keV 


NR,hc 






1 
2 


-0.62 


-1.46 


/ + 


1 
2 


_ 49 
2 


6.3 keV 


NR,hc 





9.5. Remarks on binding 

Clearly, this summary indicates that the theory of monopole binding to nuclear magnetic 
dipole moments is rather primitive. The angular momentum criteria for binding is 
straightforward; but in general (except for relativistic spin 1/2) additional interactions 
have to be inserted by hand to regulate the potential at r = 0. The results for binding 
energies clearly are very sensitive to the nature of that additional interaction. It cannot 
even be certain that binding occurs in the allowed states. In fact, however, it seems 
nearly certain that monopoles will bind to all nuclei, even, for example. Be, because 
the magnetic field in the vicinity of the monopole is so strong that the monopole will 
disrupt the nucleus and will bind to the nuclear, or even the subnuclear, constituents. 

9.6. Binding of monopole-nucleus complex to material lattice 

Now the question arises: Can the magnetic field in the detector extract the monopole 
from the nucleus that binds it? And if not, is the bound complex of nucleus and 
monopole rigidly attached to the crystalline lattice of the material? To answer the 
former question we regard it as a simple tunneling situation. The decay rate is estimated 
by the WKB formula 



Magnetic Monopoles 



59 



where the potential is crudely that due to the dipole interaction and the external 
magnetic field, 

V = -^-gBr, (9.146) 

Ai is the nuclear mass <^ monopole mass, and the inner and outer turning points, a 
and b are the zeroes of E ~ V. Provided the following equality holds, 

(9.15) 

which should be very well satisfied, since the right hand side is about 10~^^|m'|'^ MeV^, 
for the CDF field of -B = 1.5 T, we can write the decay rate as 

r ~ |m'|"^/^102^s"^exp 

where the characteristic field, defined by eBo = ml, is 4 x 10^ T. If we put in i? = 1.5 T, 
and A = 27, —E = 2.6MeV, appropriate for llA\, we have for the exponent, for 
m' = 1/2, —2 X 10^^, corresponding to a rather long time! To get a 10 yr lifetime, 
the binding energy would have to be only of the order of 1 eV. Monopoles bound with 
kilovolt or more energies will stay around forever. 

Then the issue is whether the entire Al atom-monopole complex can be extracted 
with the 1.5 T magnetic field present in CDF. The answer seems to be unequivocally 
no. The point is that the atoms are rigidly bound in a lattice, with no nearby site into 
which they can jump. A major disruption of the lattice would be required to dislodge 
the atoms, which would probably require kilovolts of energy. Some such disruption was 
made by the monopole when it came to rest and was bound in the material, but that 
disruption would be very unlikely to be in the direction of the accelerating magnetic 
field. Again, a simple Boltzmann argument shows that any effective binding slightly 
bigger than 1 eV will result in monopole trapping "forever." This argument applies 
equally well to binding of monopoles in ferromagnets. If monopoles bind strongly to 
nuclei there, they will not be extracted by 5 T fields, contrary to the arguments of 
Goto et al |154j . The corresponding limits on monopoles from ferromagnetic samples of 
Carrigan et al |155j are suspect. 

10. Searches for magnetic monopoles 

With the advent of "more unified" non-Abelian theories, classical composite monopole 
solutions were discovered, as briefiy discussed in section HI The mass of these mono- 
poles would be of the order of the relevant gauge-symmetry breaking scale, which for 
grand unified theories is of order 10^^ GeV or higher. But there are models where 
the electroweak symmetry breaking can give rise to monopoles of mass ~ 10 TeV 
jnSllMlinSlinn]- Even the latter are not yet accessible to accelerator experiments, so limits 
on heavy monopoles depend either on cosmological considerations (see for example |156j ) 
or detection of cosmologically produced (relic) monopoles impinging upon the earth or 
moon [211 ESI I157| I158| 1159^ \W{)\ \W2\ . Since the revival of interest in monopoles 
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in the 1970s, there have been two well-known announcements of their discovery: that 
of Price et al |163j . who found an cosmic ray track etched in a plastic detector, and 
that of Cabrera |158j . who reported a single event in a induction loop. The former 
interpretation was immediately refuted by Alvarez |164j . while the latter has never been 
duplicated, so is presumed spurious. 

However, a priori, there is no reason that Dirac/Schwinger monopoles or dyons of 
arbitrary mass might not exist: In this respect, it is important to set limits below the 1 
TeV scale. 

10.1. Direct searches 

In this review we will concentrate on recently obtained limits, since periodic reviews 
of the status of magnetic monopole searches have been published jl65^ I166j . Before 
2000, the best previous direct limit on magnetic monopoles was that obtained at 
Fermilab by Bertani et al jl67j who obtained cross section limits of 2 x 10~^^ cm^ 
for monopole masses below 850 GeV. As we shall see below, the Oklahoma experiment 
PT] , while not extending to as high masses, gives cross section limits some two orders of 
magnitude smaller. The recent CDF experiment ||168ij sets a three order of magnitude 
improvement over jl67j . (In contrast to |1651 I166j . we call all of these experiments 
"direct," whether they are searching for previously produced monopoles trapped in 
material or the ionization and radiation produced by monopoles passing through a 
detector.) As noted in section there have been experiments to search for monopoles 
by extracting them from matter with strong magnetic fields |155| I169j : as remarked 
there, it is doubtful such an experiment would succeed, since the binding energy of a 
monopole to the lattice is probably at least in the keV range, while the energy acquired 
by a Dirac monopole in a 100 kC field over an atomic distance is only 20 eV. 

Cosmologically produced monopoles are commonly assumed to arise from a GUT 
(Grand Unified Theory) where a grand unified group such as SU(5) breaks down into 
the standard model group SU(3) xSU(2) xU(l). Barring premature unification due to, 
say, large extra dimensions, the mass of such a monopole is expected to be of the order 
10^^ GeV, so they are incapable of being produced in accelerators. However, since in 
the early universe at least one monopole should be produced per causal domain, too 
many monopoles would have been produced |1701 1171] , and would come into conflict 
with the Parker bound, which states that cosmic fields would be quenched if the density 
of magnetic monopoles is too high |172j . This is one of the problems solved by inflation. 

Various experiments have been conducted to look for cosmic monopoles. An 
interesting limit comes from the Rubakov-Callan mechanism for monopole catalysis 
of proton decay jl73l I174j , 

M + p^M + e+ + 7r°, (10.1) 

where MACRO |175j found a limit on the flux of 3-8x10^^^ cm^^ s^^ sr^^. However, 
macro's best limit |176j . based on scintillation counters, limited streamer tubes, and 
nuclear track detectors, give a much better limit of 1.4x 10~^^ cm~^ s~^ sr~^ for monopole 
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velocities in the range 4 x 10~^ < (3 < 1, roughly a factor of two improvement over 
previous limits, and well below the Parker bound of ~ 10~^^ cm"^ s^^ sr^^. Even 
smaller limits depend on the mechanism by which a monopole would produce a track 
in ancient mica |177tll78j . One should note that lower mass monopoles, with masses of 
order 10^° GeV, arising from intermediate stages of symmetry breaking below the GUT 
scale, would not catalyze proton decay |179tll8Uj . but the more stringent MACRO limits 
still apply. 

We will discuss the three recent direct search limits in sections fTTHT!?! 
10.2. Indirect searches 

In the above, and in the following sections, we discussed direct searches, where 
the monopoles are searched for as free particles. The indirect searches that have 
been proposed and carried out rely on effects attributable to the virtual existence of 
monopoles. De Riijula in 1995 |181j proposed looking at the three-photon decay of 
the Z boson, where the process proceeds through a virtual monopole loop, as shown 
in figure El If we use his formula for the branching ratio for the Z — > 87 process, 
compared to the current experimental upper limit |182j for the branching ratio of 10~^, 
we can rule out monopole masses lower than about 400 GeV, rather than the 600 GeV 
quoted by De Riijula. Similarly, Ginzburg and Panfil in 1982 |188j and more recently 
Ginzburg and Schiller in 1999 |1841 1185j considered the production of two photons with 
high transverse momenta by the collision of two photons produced either from e~^e~ 
or qq collisions. Again the final photons are produced through a virtual monopole 
loop. Based on this theoretical scheme, an experimental limit was given by the DO 
collaboration |186j . which sets the following bounds on the monopole mass M: 



where 5* is the spin of the monopole, and m' = eg is the magnetic charge quantization 
number. 

It is worth noting that a lower mass limit of 120 GeV for a Dirac monopole has 
been set by Graf, Schafer, and Greiner |187j . based on the monopole contribution to the 
vacuum polarization correction to the muon anomalous magnetic moment. (Actually, 
we believe that the correct limit, obtained from the well-known textbook formula for 
the ^f-factor correction due to a massive Dirac particle is 60 GeV.) 

10.2.1. Difficulty with indirect limits The indirect limits mentioned above rely up the 
Feynman graph shown in figure E| If the particle in the loop is an ordinary electrically 
charged electron, this process is well-known. If, further, the photons involved are of 
very low momentum compared the the mass of the electron, then the result may be 
simply derived from the well-known Euler-Heisenherg Lagrangian |188^ 1189^ I19()j which 




for 5 = 
for S = l/2 , 
for 5 = 1 



(10.2) 
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Figure 13. 

loop. 



The light-by-light scattering graph for either an electron or a monopole 



for a spin-1/2 charged-particle loop in the presence of weak homogeneous electric and 
magnetic fields is 

^ = - T7-^" + — 7Z^[4^^')' + (10.3) 
lovr 360 (47r)^ 

where m is the mass of the particle in the loop. The Lagrangian for a spin-0 and spin-1 
charged particle in the loop is given by similar formulas, which are derived in |190t I191j 
and (imphcitly) in Ref. [US CHSl IMI , respectively. 

Given this homogeneous-field effective Lagrangian, it is a simple matter to derive 
the cross section for the 77 — *• 77 process in the low energy limit. (These results can, 
of course, be directly calculated from the corresponding one-loop Feynman graph with 
on- mass-shell photons. See |19H 1195] .) Explicit results for the differential cross section 
are given in textbooks: 

do" 139 4 / 9 / ,x 

dH = 32400^" V^ + "^*>' <^°-*' 
and the total cross section for a spin-1/2 charged particle in the loop is 

973 

^=T7^T^"'— ' ^/m<^l, s = Au;^ (10.5) 
101257r m* 

The numerical coefficient in the total cross section are 0.00187, 0.0306, and 3.50 for spin 
0, spin 1/2, and spin 1 particles in the loop, respectively. 

How is this applicable to photon scattering through a monopole loop? This would 
seem impossible because of the existence of the string, which renders perturbation theory 
meaningless. Of course, no one has attempted a calculation of the "box" diagram with 
the monopole interaction. Rather, De Riijula and Ginzburg (explicitly or implicitly) 
appeal to duality, that is, the dual symmetry ()2.2a|) that the introduction of magnetic 
charge brings to Maxwell's equations: 

E^B, B^-E, (10.6) 

and similarly for charges and currents. Thus the argument is that for low energy photon 
processes it suffices to compute the fermion loop graph in the presence of zero-energy 
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photons, that is, in the presence of static, constant fields. Since the Euler-Heisenberg 
Lagrangian is invariant under the duahty substitution on the fields alone, this means we 
obtain the low energy cross section a^^^^^ through the monopole loop from the equation 
for the QED cross section by the substitution e ^ g, or 

a^ag = 137m'^, 2|m'| = 1, 2, 3, . . . . (10.7) 

It is critical to emphasize that the Euler-Heisenberg Lagrangian is an effective 
Lagrangian for calculations at the one fermion loop level for low energy, i.e., u/m ^ 1. 
However, it becomes unreliable if radiative corrections are large. (The same has been 
noted in another context by Bordag, Robaschik, Wieczorek, and Lindig |196tll97j .) For 
example, the internal radiative correction to the box diagram have been computed by 
Ritus HlHI and by Renter, Schmidt, and Schubert MI in QED. In the 0{a^) 
term in the expansion of the EH Lagrangian (fTTOIl . the coefficients of the (F^)^ and the 
[FFy terms are multiplied by 

4:0a ^, o,\ , / 1315a ^, o,\ 
l + y- + 0(a')j and (^1 + -_- + 0(^2) j , (10.8) 

respectively. The corrections become meaningless when we replace a ^ ag. 

10.2.2. Unitarity hound This would seem to be a devastating objection to the results 
given by Ginzburg et al |184| 1185] and used in the DO analysis |186j . But even if one 
closes one's eyes to higher order effects, it seems clear that the mass limits quoted are 
inconsistent. 

If we take the cross section given by ()10.5|) and make the duality substitution, we 
obtain for the low energy light-by-light scattering cross section in the presence of a 
monopole loop (M is the monopole mass) 

If the cross section were dominated by a single partial wave of angular momentum J, 
the cross section would be bounded by 

,<^P£±I).!^, j.i. (10.10) 

s s 

Comparing this with the cross section given above, we obtain the following inequality 

for the cross section to be consistent with unitarity, 

M , „ , , 

— >6|m'|. (10.11) 

UJ 

But the limits quoted by DO for the monopole mass are less than this: 
M 



, - > 870 GeV, spin 1/2, (10.12) 
2\m'\ 

because, at best, a minimum estimate is {oj) ~ 300 GeV, so the theory cannot sensibly 
be applied below a monopole mass of about 1 TeV. (Note that changing the value of 
J in the unitarity limits has very little effect on the bound since an 8th root is taken: 
Replacing J by 50 reduces the limit only by 50%.) 
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Similar remarks can be directed toward the De Riijula limits |181j . That author, 
however, notes the "perilous use of a perturbative expansion in g." However, although he 
writes down the correct vertex, he does not, in fact, use it, instead appealing to duahty, 
and even so he admittedly omits enormous radiative corrections of 0{ag) without any 
justification other than what we believe is a specious reference to the use of effective 
Lagrangian techniques for these processes. 

As we will see, some of the same objections apply to the direct search limits. 
The advantage, however, of the latter, is that the signal of a positive event is more 
unambiguous, and in the Oklahoma and HI experiments, a monopole, if found, would 
be available for further study. 



11. Oklahoma experiment: Fermilab E882 

The best prior experimental limit on the direct accelerator production of magnetic 
monopoles is that of Bertani et al in 1990 jl67j (see also Price et al |20H 1202] ): 

(T < 2 X IQ-^'^cm^ for a monopole mass M < 850 GeV. (11.1) 

The fundamental mechanism is supposed to be a Drell-Yan process, 

p + p ^ M + M + X, (11.2) 



where the cross section is given by 
- (68. Sn^)^/?"* 



dM ' ' ' 9s 



f^j:Q'M^^)^^(^). (11.3) 

J Xi \SXi J 



Here Ai is the invariant mass of the monopole-antimonopole pair, and we have included 
a factor of to reflect (1) phase space and (2) the velocity suppression of the magnetic 
coupling, as roughly implied by (j5.5|) -See also (j2.16|) . Note that we are unable to 
calculate the elementary process 

qq^l* ^ MM 

perturbatively, so we must use nonperturbative estimates. 

Any monopole produced at Fermilab is trapped in the detector elements with 100% 
probability due to interaction with the magnetic moments of the nuclei, based on the 
theory described in section IHl The experiment consists of running samples obtained 
from the old DO and CDF detectors through a superconducting induction detector. 
Figure El is a sketch of the DO detector. We are able to set much better limits than 
Bertani et al |167j because the integrated luminosity is 10"^ times that of the previous 
1990 experiment: 

/: = 172±8pb"^ (DO). (11.4) 

We use energy loss formula of Ahlen |133| I132| |203|, as described in section |H1 The 
graph in figure IT^ shows the energy loss dE/dx for various materials. 
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Chamber QiambCT Railialion ChainlKr 
Deuctor 

Fig. 3. ArraRgcment of the DO traddng and tiaiisilian radiatiofi detecton. 



Figure 14. Arrangement of the DO tracking and transition radiation detectors. 



Figure ^1 is a diagram of the OU magnetic monopole induction detector. It is a 
cylindrical detector, with a warm bore of diameter 10 cm, surrounded by a cylindrical 
liquid N2 dewer, which insulated a liquid He dewer. The superconducting loop detectors 
were within the latter, concentric with the warm bore. Any current established in the 
loops was detected by a SQUID. The entire system was mechanically isolated from the 
building, and magnetically isolated by fi metal and superconducting lead shields. The 
magnetic field within the bore was reduced with the help of Helmholtz coils to about 
1% of the earth's field. Samples were pulled vertically through the warm bore with a 
computer controlled stepper motor. Each traversal took about 50 s; every sample run 
consisted of some 20 up and down traversals. Most samples were run more than once, 
and more than 660 samples of Be, Pb, and Al from both the old CDF and DO detectors 
were analyzed over a period of 7 years. 



11.1. Monopole induced signal 

Note that if the shield were not present, the supercurrent induced by a monopole of 
strength g passing through a loop of radius r and inductance L would be given by 

where z{t) is the vertical position of the monopole relative to the position of the center 
of the loop. A more detailed theory is described in the following. The theory can be 
verified with a pseudopole, which is a long, ~ 1 m, electromagnetic solenoid, which 
produces a field near one end very similar to that produced by a pure magnetic pole. 
The excellent agreement between theory and experiment is indicated in figure E| 
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Figure 15. Sketch of the OU induction detector. Shown is a vertical cross section; 
it should be imagined as rotated about the vertical axis labelled "centerline." 
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Figure 16. Typical step plots: DO aluminum, CDF lead, and CDF aluminum. The 
experimental data was collected from pseudopole simulations; the steps shown are 
for the difference between the results with reversed polarizations of the pseudopole. 
Data agrees well with the theory which incorporates the effect of the shielded 
superconducting loops. The theory without the shield, given by Barger and OUson 
|2l)4| . is also shown. 
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Figure 17. Diagram of monopole detector. Tiie monopole g is assumed to be on 
the central axis at a height z' above the bottom of the detector, which we model 
as a cylindrical perfectly conducting can of radius a, closed at the bottom. The 
superconducting loop of radius r is a height Z above the base. 



11.1.1. Simplified theory of monopole detector This subsection describes the basis of 
the functioning of our magnetic monopole detector. It works by detecting the magnetic 
flux intercepted by a superconducting loop contained within a superconducting cylinder. 
The detector is sketched in figure IT71 

In order to incorporate finite-size effects, we consider first a perfectly conducting 
right circular cylinder of radius a of semi-infinite length, with axis along the 2;-axis, and 
with a perfectly conducting circular bottom cap at 2; = 0. We use cylindrical coordinates 
p, 6, and z. 

Because the boundaries are superconductors, the normal component of B must 
vanish on the surfaces, that is, 



= 0, 

p — a 

2>0 



= 0. (11.6) 

=0 



Now suppose a magnetic pole of strength g is placed on the z axis at z = z' > 0. 
This could either be a magnetic monopole (magnetic charge) or one pole of a very 
long electromagnet ("pseudopole"). Imagine a circular conducting loop of radius r < a 
centered on the axis of the cylinder and perpendicular to that axis, with center at 
z = Z. Inside the cylinder and outside of the loop, B is derivable from a magnetic 
scalar potential, 

B = -V0M, (11.7) 

since we may ignore the displacement current, because the time variation is negligible. 
(pM satisfies Poisson's equation, in cylindrical coordinates: 
/Id d 1 92 52 \ 
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where r' is the position of the monopole, r' = {p',6',z'). This is the equation for a 
Green's function, which we can express in separated variables form. That is, we write 

2 



TT 



/ dkcoskzcoskz' y — e^"^(^-^')^„(p,p';A:), (11.9) 
Jo _ 27r 



where, in view of the first boundary condition in (jll.6p . we may express the reduced 
Green's function in terms of modified Bessel functions: 



gm{p,p'] k) = -ATTgI^{kp, 



;ii.io) 



where p< (p>) is the lesser (greater) of p, p'. If the monopole is confined to the z axis, 
only the m = term survives: 



— - f dk cos kz cos kz' 
Jo 



Koikp) + loikp) 



Ki{ka) 
Ii{ka) 



which uses 

J^(x) = Ji(x), K',{x) = -K,{x). 
By integrating over the cross section of the loop using 



dttKoit) 



-X KAx) + 1, 



dtt Io(t) = X Ii{x) 



we obtain the following formula for the magnetic flux subtended by the loop, 

jdS-B = 47rg [r]{Z - z') ~ F{Z, z')] , 

where the step function is fl3.31aj) . and the response function is 

K,{x) 



(x) 



(11.11) 
(11.12) 
(11.13) 

(11.14) 
(11.15) 



2r f°° z z' \ K] 

F(z,z') = / dx sinx- cosx— < Ki(xr/a) — Ii(xr /a) — 
7!' a Jq a a I Ii 

Now suppose that the pole is slowly moved from a point far above the loop, 
z' = +00, to a point below the loop, z' = zq, Z > zq. Then from Maxwell's equation 



VxE = — — B J^, 

cot c 

where Jm is the magnetic current density, the emf induced in the loop is 

d<l> Att 



E-dl,= -— + — 55(t), 
cdt c 



:ii.i6) 



:ii.i7) 



if t = is the time at which the pole passes through the plane of the loop. The net 
change in emf gives rise to a persistent current / in the superconducting loop. 



LI 



1 At: An 

£dt = — A$ + —g = —gF{Z, zq), 
c c c 



:ii.i8) 



where L is the inductance of the loop, and the response function F is given in (|11.15|) . 
This is just a statement of the Meissner effect, that the flux change caused by the moving 
monopole is cancelled by that due to the current set up in the loop. 



Magnetic Monopoles 



69 



When the loop is very far from the bottom cap, Z S> a, only small x contributes 
to the integral in Eq. ()11.15|) . and it is easy to see that 

£dt = ^ (l (11.19) 
c V « / 

so the signal is maximized by making the loop as small as possible, relative to the 
radius of the cylinder. We get the full flux of the monopole only for a loop in empty 
space, a/r — > oo. This perhaps counterintuitive effect is due to the fact that the 
superconducting walls confine the magnetic fiux to the interior of the cylinder. Thus 
for the superconducting can, the induced current in the detection loop caused by the 
passage of a monopole from 2;' = oo to 2;' = is 

Atto A$ Atto ^(z' = 0) , , 

LI = — = — ^ ^, 11.20 

c c c c 

which yields the result flll.l9j) if one assumes that the magnetic field is uniform across 
the can's cross section at the position of the loop when the pole is at the bottom, because 
all the fiux must pass up through the can. If we consider, instead, an infinite, open- 
ended, superconducting cylinder, with the monopole passing from z = +00 to z = —00, 
at either extreme half the fiux must cross the plane of the loop, so with the uniformity 
assumption we get the same result: 

L/ = ^-^ = ^ri-Cl. m.2ii 



c c c \ 



The simple assumption of a uniform magnetic field is apparently justified by the exact 
result (|11.19|1 . 

We conclude this discussion by noting how the exact calculation is modified for an 
infinite superconducting cylinder. In the magnetic scalar potential, the integral over k 
mode functions in ()11.9|) is replaced by 

"^^-e'K^-^'l (11.22) 



27r 

which has the effect of replacing the flux expression (|11.14p by 

$ = 2TTg[e{z - z') - F{Z - z', 0)], (11.23) 

where e(^) is given by ()3.315|) . Then the induced current in the detection loop when the 
monopole passes from a point above the loop z' = Z + ^ to a. point, equidistant, below 
the loop, z' = Z — ^, is 

L/ = i^f«.0)^l=«(l-4). (11.24) 

c c \ a^/ 

where the last limit applies if ^/a ^ 1. This result coincides with that in (jll.l9|) . The 
function R{^) = 0)/(l — r"^ ja?) + i, corresponding to a monopole starting from a 

point Zi far above the loop, zi — Z ^ a, and ending at a point zq = Z — S,-, is plotted 
as a function of ^ for our parameter values in flgure IT^ . where it is shown to agree well 
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with experimental data. This response function coincides with the result obtained from 
(jll.lSp . because 

FiZ, Z-0 = \f{2Z - 0) + \F{i, 0) ^ 1 (^1 - + 0), (11.25) 

ii Z/a ^ 1. This shows that the effect of the endcap (which of course is not present 
in actual detector) is negligible, demonstrating that the fact that the superconducting 
shield is of finite length is of no significance. 



11.2. Background effects 

All nonmagnetic but conducting samples possess: 

• Permanent magnetic dipole moments /x, which give rise to signals in free space of 
the form 



m 



[11.26) 



Lc [r2 + ^(t)2]3/2 

Induced magnetization: Conducting samples passing though magnetic gradients 
with speed v produce time-varying magnetic fields which induce signals in our 
detector, 

N V 1 f , , . o , dBz , ,.1 f z' a 
nt) = -J / idryair')—^{z')-H -, 

C'^ L J oz r \r r ^ 

where H is the response function, essentially that appearing in (jll.lSp . 



;il.27) 



dy y cosy— 



Ki{ya/r) 
h{ya/r) _ 



2 (r2 + ;2'2^3/2' 



a/r 



oo. 



;il.28a) 
(11.286) 



11.3. Calibration, real data, and limits 



The pseudopole data shown in figure UHl clearly shows that we could detect a Dirac 
pole. We demonstrated that the detector (SQUID response) was remarkably linear over 
a range of 0.7-70 Dirac poles. 

As one sees from figure ^1 real samples have large dipole signals; what we are 
looking for is an asymptotic step indicating the presence of a magnetic charge. Steps 
seen are typically much smaller than that expected of a magnetic pole of Dirac strength. 
The histograms of steps are shown in figures IT^HTTl 

For m' = 1/2 the 90% confidence upper limit is 4.2 signal events for 8 events 
observed when 10 were expected |205j . These 8 samples were remeasured and all fell 
within ±1.47 mV of m' = 0. (More than 1.28o- from \m'\ = 1/2.) For m' = 1 the 90% 
confidence upper limit is 2.4 signal events for zero events observed and zero expected. 

By putting in angular and mass acceptances we can get cross section limits as shown 
in table El These numbers reflect the new analysis, published in 2004 |2I], and so differ 
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Figure 18. Steps: DO Al, CDF Pb, and CDF Al. 
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Figure 19. Steps from DO samples. A Dirac pole would appear as a step at 2.46 mV. 

somewhat from our earlier published results [20] . To obtain the mass limits, we use the 
model cross sections given in figure |221 

Finally, we show in figure ESI what might be achievable at the LHC, using the same 
techniques applied here. 
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Figure 20. Steps from CDF Pb samples. A Dirac pole would appear as a step at 2.46 
mV. 
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Figure 21. Steps from CDF Al samples. A Schwinger pole {^gi)) would appear as a 
step at 10.64 mV. 



12. HI limits 

We now turn to the limits on monopole production obtained from e^p collisions at 
HERA recently published by the HI collaboration |2U6j . This production mechanism 
is intermediate between that of pp experiments such as that given in |2I], and that of 
the possible production through e~^e~ collisions |182j . and conceivably might yield a 
cleaner interpretation if monopole condensates are responsible for the confinement of 
quarks ^SJ I2U71 12U81 l2U9j . Although the mass limits determined are not as strong as 
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Table 2. Alternative interpretations for different production angular distributions of 
the monopoles, comparing 1 and f ± cos^ 9. Here the cross section aa corresponds to 
the distribution 1 + a cos^ 6, and similarly for the mass limits (all at 90% confidence 
level) . 
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Figure 22. Cross section vs. mass limits. The three graphs show three different 
assumptions about the angular distribution, since even if we knew the spin of the 
monopole, we cannot at present predict the differential cross section. Shown in the 
second figure are the Bertani |167| and lunar limits. 
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Figure 23. Monopole pair masses as a function of the cross section at the Tevatron 
{pp at 2 TeV) and at the LHC {pp at 14 TeV). Both include the (3'^ correction, and are 
for a Dirac monopole, m' — 1/2. 



in our experiment described in the previous section, it is crucial that different physical 
domains be explored carefully. 

In their experiment, the aluminum beam pipe used at the HI interaction point at 
HERA during 1995-1997 was cut into 75 long and short strips. This beam pipe had 
been exposed to an integrated luminosity of 62 ± 1 pb~^. These strips were then placed 
on a conveyor belt and passed through a warm-bore magnetometer at Southhampton 
Oceanographic Centre, UK. If a monopole passed through the superconducting coil, as 
in our experiment, it would establish a persistent current there, which would be detected 
by a SQUID. Again they calibrated their detector by a long, thin solenoid, which at 
each end produced a pseudopole. Calibration was within some 10% for pole strengths 
above go = hc/2e. Large dipole signals were seen, but the signals always returned to 
the baseline unless a pseudopole was present. A few persistent current events were seen, 
but they always disappeared upon remeasurement. Some of the runs exhibited large 
fluctuations of unknown orgin, but none was consistent with a monopole event. 

No monopole was detected in their experiment of strength greater than O.lgo for 
a sample consisting of 93 ± 3% of the beam pipe. To interpret this as a limit on the 
production cross section, models had to be adopted, since perturbation theory was 
unreliable. Two models were tried: 



where in both models the monopole pairs were produced by two-photon processes. 
The effects on the produced monopoles by the HI magnetic fields were included. The 
stopping power was computed by using the classical results of Ahlen |210| I203| I132j . 

The results of their analysis are expressed in plots of the upper limits on the cross 
sections for a given monopole mass, up to a mass of 140 GeV, for different models and 



e+p e^MMp, 
e+j9 ^ e+MMX, 



spin monopole, 
spin 1/2 monopole. 
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magnetic charges. That is, cross sections above those hmits are excluded based on their 
experimental analysis. These limits are weakest for the Dirac charge, g^, and strongest 
for Schwinger quantization based on quark charges, Qgo, of course. For model ()12.1aj) 
the cross section limits range from about 1 pb for several GeV to more that 100 pb for 
140 GeV for g^- For higher charges the limits are relatively constant around 0.1 pb or 
less. For ()12.16|) the limits are similar, except for gjy, where the limit drops below 0.1 
pb for 10 GeV masses or less. Thus their results our complementary to ours: Our mass 
limits are stronger, but in some cases they reach smaller cross sections for lower masses. 

13. CDF limits 

Quite recently, a new CDF experiment |168j has been announced, which claims, on 
the basis of an integrated luminosity of 35.7 pb~^, a production cross section hmit for 
spin-1/2 monopoles below 0.2 pb for masses between 200 and 700 GeV, and hence, in 
a Drell-Yan model, a lower mass limit of 360 GeV. (These limits are quoted at the 95% 
confidence level.) This is based on quite different technology than the Oklahoma or HI 
experiments. Rather, they looked at a sample of pp events collected during 2003 by the 
CDF detector by a special trigger. The signal for a monopole is the large ionization 
and heavy production of delta rays by such a particle. They use our crude model of 
replacing e by g(3 in the Drell-Yan production mechanism [2011^; apart from this simply 
replacing the lepton mass by the monopole mass. Acceptance is affected by production 
kinematics, which effect they estimate at 10%. Light monopoles will be swept out of 
the detector by the magnetic field, while heavy monopoles may reach the time-of-fiight 
detector too late to cause a trigger. Other particles ("spoilers") may cause a charge 
integration to start in the detector before a monopole signal arrives; they estimate a 
few percent fraction for a monopole of mass 400 GeV. 

Out of 130,000 candidate events, no monopole trigger events were found, from which 
the limit quoted above was extracted. They believe they can push the limit on masses 
up another 100 GeV with additional Run II data. 

13.1. Comments on CDF experiment 

One might ask how can CDF claim stronger limits than we do based on less than 
1/4 of the integrated luminosity of our experiment. The answer, I believe, is that 
our experimental limit is extremely conservative. They may have underestimated the 
systematic effect of huge uncertainties in the production mechanism, while at the 
same time they claim our limit is dependent on the trapping model, which it is not. 
Undoubtedly, the dE/dx signature of monopoles is much less well understood than the 
clear-cut electromagnetic signature of an induction detector. This is not to denigrate 
the utility of this measurement, but to emphasize that the limits so obtained are subject 
to large, relatively uncontrolled, uncertainties. 
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14. Conclusions 



One magnetic monopole, carrying magnetic charge g, will result in the quantization of 
electric charge throughout the universe, 
m'hc 



9 

where m' is a half-integer, 



(14.1) 



m' = 0,±^,±l,±^,±2,±^,.... (14.2) 

That electric charge is quantized in integer multiples of the electron charge (or integer 
multiples of the quark charge) is an overwhelming fact, which does not possess a simple 
explanation. This, perhaps, is the most compelling argument in favor of the existence 
of magnetic charge. A second argument is the greater symmetry (duality) imparted to 
Maxwell's equations, and to classical and quantum electrodynamics, if both electric and 
magnetic charges are present. Thus from phenomenological and theoretical bases, the 
arguments in favor of the existence of magnetic charge and for dual QED are at least as 
strong as those for supersymmetry. Unfortunately, like for the latter, there is not a real 
shred of observational evidence in favor of magnetic charge, although here it is far less 
embarrassing to be in that situation, since the most likely mass range for magnetically 
charged particles is not far from the Planck scale. 

In this review we have concentrated on the theory of point Dirac monopoles or 
Schwinger dyons, starting from the classical scattering, through the nonrelativistic 
quantum mechanical description, to the quantum field theory of such objects. For 
lack of space we have only briefly referred to the classical monopoles that arise from the 
solution of non-Abelian gauge theories. From the point of view of phenomenology and 
the setting of experimental limits, the point description should be adequate, since the 
structure of composite monopoles only emerges at the energy scale that sets the mass of 
the particles. (An exception, of course, occurs with limits based on that structure, such 
as the catalysis of proton decay.) In addition, our concern has been chiefly with the 
quantum description, which has been only roughly sketched for composite monopoles. 

We close, as did Schwinger in his provocative article [llj, by quoting from Faraday: 
"Nothing is too wonderful to be true, if it be consistent with the laws of nature, and in 
such things as these, experiment is the best test of such consistency." 
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